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Abstract 

Dimensional regularization is applied to the computation of the gravitational wave field generated 
by compact binaries at the third post-Newtonian (3PN) approximation. We generalize the wave 
generation formalism from isolated post-Newtonian matter systems to d spatial dimensions, and 
apply it to point masses (without spins), modelled by delta-function singularities. We find that 
the quadrupole moment of point-particle binaries in harmonic coordinates contains a pole when 
e = d — 3—>0at the 3PN order. It is proved that the pole can be renormalized away by means 
of the same shifts of the particle world-lines as in our recent derivation of the 3PN equations of 
motion. The resulting renormalized (finite when e ^ 0) quadrupole moment leads to unique values 
for the ambiguity parameters C ^ C which were introduced in previous computations using 
Hadamard’s regularization. Several checks of these values are presented. These results complete 
the derivation of the gravitational waves emitted by inspiralling compact binaries up to the 3.5PN 
level of accuracy which is needed for detection and analysis of the signals in the gravitational-wave 
antennas LIGO/VIRGO and LISA. 
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I. INTRODUCTION 


A compelling motivation for accurate computations of the gravitational radiation field 
generated by compact binary systems {i.e., made of neutron stars and/or black holes) is 
the need for accurate templates to be used in the data analysis of the current and future 
generations of laser interferometric gravitational wave detectors. It is indeed recognized 
that the inspiral phase of the coalescence of two compact objects represents an extremely 
important source for the ground-based detectors LIGO/VIRGO, provided that their total 
mass does not exceed say 10 or 20 Mq (this includes the interesting case of double neutron- 
star systems), and for space-based detectors like LISA, in the case of the coalescence of two 
galactic black holes, if the masses are within the range between say 10^ and 10® M©. 

For these sources the post-Newtonian (PN) approximation scheme has proved to be the 
appropriate theoretical tool in order to construct the necessary templates. A program started 
long ago with the goal of obtaining these templates with 3PN and even 3.5PN accuracy.^ 
Several studies have shown that such a high PN precision is 

probably sufficient, not only for detecting the signals in LIGO/VIRGO, but also for analyzing 
them and accurately measuring the parameters of the binary (such high-accuracy templates 
will also be of great value for detecting massive black-hole mergers in LISA). The templates 
have been hrst completed through 2.5PN order, for both the phase |llLll2Lll3L ll4| and wave 
amplitude ElEl. The 3.5PN accuracy for the templates (in the case where the compact 
objects have negligible intrinsic spins) has been achieved more recently, in essentially two 
steps. 


(1) The hrst step has been to compute all the terms 


motion, either in Hamiltonian form |17L [1 


Q El in 0, and the 3.5PN 


generation formalism SEIE 


in both the 3PN equations of 
or using harmonic coordinates 
yavitational radiation held, using a multipolar wave 
l28j |. by means of the Hadamard self-held regulariza¬ 
tion j29|, |30|, |3l|, [sl 5 ill short HR. (The 3.5PN terms in the equations of motion have 
been added in Refs, m El; El-) However, a few terms were left undetermined by 
Hadamard’s regularization, which corresponds to some incompleteness of this regu¬ 
larization occurring at the 3PN order. These terms could be parametrized by some 
unknown numerical coefficients called ambiguity parameters. 

(2) The second step has been to hx the values of the ambiguity parameters by means of 


dimensional regularization j36|, |37|, |38[ , henceforth abbreviated as DR. Technically, DR 


is based on analytic continuation in the dimension of space d = 3 -f £. The ar nbiguity 
parameter A entering the 3PN equations of motion has been computed in Refs. (391 l4n| | , 
with result A = —1987/3080. (This result has also been obtained with an alternative 
approach in Refs. (dll . Idl . Idl] .) The three ambiguity parameters appearing in the 3PN 
gravitational radiation held will be shown in the present paper to have the following 
unique values 

9871 7 

e = -^, ft: = 0, ( = (1.1) 

s g 24 g ) , s 22 , V y 

as already announced in Ref. g. The method we use for applying DR essentially 
consists in computing the difference between DR and some appropriately dehned 


^ Following the standard custom we use the qualifier nPN for a term in the wave form or (for instance) the 
energy flux which is of the order of 1/c^" relatively to the lowest-order Newtonian quadrupolar radiation. 


2 











Hadamard-type regularization called below the pure-Hadamard-Schwartz (pHS) reg¬ 
ularization. 

Those results complete the determination of the 3.5PN-accurate phase evolution as it suffices 
to insert into the formulas of Ref. the value for A, together with the values given by 
(mi). Actually, this phase evolution depends only on A and on the following particular 
combination of parameters, 


9 = f + 2K + C = -. (1.2) 

The present paper is devoted to the details of our DR computation of the ambiguity param¬ 
eters, item (2) above, which has led to the values (j1 ■2fl . We refer to for a summary 

of our method and a general discussion. 

Let us emphasize that the values (CH), which constitute the end result of the application 
of DR, have all been conhrmed by alternative methods. Our first independent check has 
been the conhrmation of one particular combination of the ambiguity parameters, namely 
^ + K, which was shown to follow from the requirement that the 3PN mass dipole moment 
of the binary, computed in from the multipolar wave generation formalism, should 
agree with the 3PN center-of-mass position, known from the conservative part of the 3PN 
equations of motion in harmonic coordinates ji^. Secondly, we have also obtained the 
value of C by considering the limiting physical situation of a boosted Schwarzschild solution, 
corresponding to the case where the mass of one of the particles is exactly zero, and the other 
particle moves with uniform velocity |4^. It can be argued from this calculation that the 
value of C in Eq. dn} is a consequence of the global Poincare invariance of the multipolar 
wave generation formalism. Thirdly, in Sec. IVIII below, we shall be able to show that the 
value of K is zero by a diagrammatic approach (where the “diagrams” are taken in the sense 
of li^), showing that no dangerously divergent diagrams contributing to k, appear at this 
order. Those checks altogether provide a conhrmation, independent from DR, for all the 
parameters (HH). 

The plan of this paper is as follows. In Section O we investigate the symmetric-trace- 
free (STF) multipole decomposition in d dimensions for a scalar held with compact-support 
source. In Section nm we generalize to d dimensions the known results for the multipole 
expansion of the gravitational held and the dehnition of the source-type multipole moments. 
Section HU is devoted to the explicit expressions of the source terms in the latter source 
multipole moments at the 3PN order in terms of a convenient set of retarded-like elementary 
potentials. Then, in Section m we obtain a general formula for the diherence between DR 
and HR (in the pHS variant of it). This diherence is non-zero at the 3PN order because of 
the occurrence of poles in d dimensions (he., oc 1/e). In Section|VTlwe deduce the ambiguity 
parameters from the DR regularization of the 3PN mass quadrupole moment, and we check 
that the 3PN mass dipole is in agreement with the known center-of-mass position deduced 
from the equations of motion. Section IVlfl deals with a direct computation of the pole part 
of the moments using diagrams, their renormalization using shifts of the world-lines, and the 
check that k = 0. In Section IVlIIl we present an alternative derivation of the value of ( based 
on considering the physical situation of a single boosted point particle in d dimensions (the 
result agrees with the recent computation of the boosted Schwarzschild solution in |45l|i. 
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II. MULTIPOLE EXPANSION OF A SCALAR FIELD IN d DIMENSIONS 


A crucial input for the derivations we are going to perform in the present article is the 
multipolar expansion of solutions of flat space-time wave equations in D = d + 1 dimensions. 
We denote by □ = the flat d’Alembertian operator, using the signature “mostly 

plus”, he., □ = A — c~‘^d^, where dt = d/dt and A is the Laplace operator. We hrst consider 
the case of a scalar wave equation, say 

□(p(x,f) = ^(x,f), (2.1) 

and shall postpone to Sec. Illll the case of tensorial wave equations. Note that, in the present 
work, we shall not introduce any numerical factor in the “source” S on the right hand side 
(RHS) of the inhomogeneous scalar wave equation (12.1|) . Similarly, we define the scalar 
Green functions as the solutions of 


□G(x,f) = (2.2) 

where (5*''^^(x) is a d-dimensional Dirac distribution, such that J f{x) = /(O). 

When d = 3, the retarded Green function takes the simple form 


^(3+1) 

'-^Ret 


(x,f) 


|x|/c) 
Att IxI 


(2.3) 


Because of the presence of the factor — l/dvr in (I2.d|l . it was convenient, when working in 
3-1-1 dimensions, to introduce a factor —dvr in front of the RHS’s of (EH) and (E2D. However, 
there is no analogous, universally simplifying factor in D dimensions, so it is finally simpler 
to introduce no factors at all in (EIII)-(E3)- 

The D-dimensional retarded Green function has no simple expression in (f, x) space. 
However, starting from its well-known Fourier-space expression, one can write the following 
simple integral expression (see e.g. jd^L 


GRet(x,f) 


Ojt) 

(27r)'^/2 




sm{ckt) Jd_^{kr) . 


(2.4) 


Notice that this is in fact a function of t and r = |x| only: say GRet(x, f) = 

Here 6 {t) is the Heaviside step function, and Jd_^(fcr) the usual Bessel function. Actually, 
we shall never need to use the explicit form (Ell of the Green function in D dimensions. 
Indeed, we shall obtain the d-dimensional generalizations of the 3-dimensional relativistic 
multipole moments, obtained in Refs, j^, , by working directly with the source S of 

the wave equation EH, or of its tensor generalizations. To do this, we note first that the 
retarded solution of m reads 


V7(x,f) 


/ 


(i‘‘y ds GRe,(x - y, i - s) S(y, s). 


(2.5) 


In this section, we shall consider sources S'(x, f) having a spatially compact support in 
d space dimensions: say S{'x,t) = 0 when |x| > a, where a is the source’s radius. We are 
interested in the multipolar expansion of the held </9(x, f), he., its decomposition (when con¬ 
sidered in the external domain |x| > a) in d-dimensional spherical harmonics. Traditionally, 
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the multipolar expansion of (p(x, t), Eq. ()2.5|1 . is obtained by expanding the spatial kernel 
GRet(x — y) in powers of |y| —0. This introduces the (reducible) multipole moments of 
the source, say f d'^yy'’^ ■ ■ ■y''^ S{y). A simpler, formally equivalent way of proceeding is to 
replace the continuous source *S'(x) by its ^^distributional skeleton’’', he., an expansion in in¬ 
creasing derivatives of the d-dimensional Dirac distribution 5(x). [For notational simplicity, 
we henceforth suppress the superscript (d) on (5(x).] This skeletonized version of the source 
S is equivalent to a continuous function S'(x) with compact support when (and only when) 
it is integrated by a regular kernel A'(x, y), as in (12.dll . It reads 

+°° (—Y 

S=“(x,() = 5^i^Si(«)SL5(x), (2,6) 

£=0 

where the coefficients are the reducible multipole moments 

Slit) = j d^yyLS{y,t). (2.7) 

We recall our simplihed notation: L denotes a multi-index i^ - ■ -ii and we use the shorthands 
• di^, where di = dldx\ and yL = yi^--- yi^, where yi = y\ 

The skeleton expansion (ESI) does not yet give rise to a multipole expansion because 
the various terms on the RHS of (ED) do not correspond to irreducible representations 
of the d-dimensional rotation group 0{d). However, it is relatively simple to transform 
the expansion (Ei into irreducible components. To do this, it is enough to decompose 
the symmetric tensors Sl into irreducible symmetric and trace-free (STF) pieces, which is 
easily done by using the STF decomposition of in d dimensions, obtained by recursively 
separating the traces, like in ytj = ^ djj |yp. Here we denote the STF projection by 

means of a hat: tjl = STF[|/ij ■ ■ -yi^], or sometimes by means of brackets surrounding the 
indices: yi = y{L)- The general formula dehned by this recursion has already been given in 
Ref. and reads 


yi 

with a\ 


"{n*2''' ®*2fe-i*2fe yL—2K} |y I ) 

A:=0 


1 r(f+ £-2fc) 

¥ r{l + e-k) 


(2.8a) 


(2.8b) 


Here, 6 ij is the Kronecker symbol, [|] denotes the integer part of |, L — 2K is a multi-index 
with i — 2k indices, and T is the usual Eulerian function. The curly brackets surrounding 
the indices refer to the (unnormalized, minimal) sum of the permutations of the indices 
which keep the object fully symmetric in L, for instance d{pl4} = dijVk + SikVj + SjkVi (for 
convenience we do not normalize the latter sum). 

We replace the STF decomposition (El) into (El and insert the resulting moments back 


We refer to the Appendix B of j4nl | for a compendium of formulae for working in a space with d dimensions. 
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(2.9a) 


into Eq. (ESD- After some simple manipulations we arrive at 


where 


+ 00 / +00 

e=o ■ k=o 

1 r(f + £) 
2^>^k\T (l + i + k) ■ 



A^L|yr"^(y,t) 


(2.9b) 


At this point let us notice that any term in the skeletonized source S'®'^®*(x, t) which is in the 
form of a d’Alembert operator □ acting on spatial gradients or time derivatives of the delta 
function, say □[95(x)],^ will give no contribution to the multipole expansion of (p(x,f). 
Indeed, a term in the source of the form [fit) 52,(5(x)], with i > 0, i > 0, will yield a 
contribution to the solution of the form (□*+^ [f{t) ^^^(x)]) = [f{t) 9l5(x)]. Such a 

contribution is localized at the spatial origin x = 0 and thus vanishes outside of the world 
tube r < a containing the source. 

We now transform the Laplacians in (E3D into d’Alembertians using 




□ + 




A k\ 





( 2 . 10 ) 


We then arrive at an irreducible (STF) decomposition of the skeletonized source S, which 
is of the type 

+00 ,Y 

5Skel(x,t) = ^Lit) dLSi^) + 0(086). (2.11) 

e=o 

Here the last term, symbolically denoted O {Od 5), is an (infinite) sum of terms of the form 
□*’''^[/(f)9i5(x)] with z > 0, £ > 0. As we just said, these terms will not contribute to the 
multipole expansion of the held 93(x, f), he., considered in the external domain r > a. 

The most useful result for our purpose is the explicit expression of the STF moments in 
Eq. ()2.11|1 which we hnd to be 


Mt) 


= / d^yyLSe(y,t), 


( 2 . 12 ) 


where we have introduced a convenient ^-dependent weighted time average given by the 
formal inhnite PN series 



(2.13) 


The coefficients af are those which have been introduced in Eq. (I2.9bjl . When written out 
explicitly, the “effective” source Si(y, t) reads. 


Si(y,t) 


1 

{2k)\\{2i + d){2i + d + 2)---(2i + d + 2k-2) 


(2.14) 



2 k 

S(y,t) + --- 


^ Here the notation d symbolizes any product of space or time derivatives (so that, for instance, d can 
involve any power of the box operator □ itself) 
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where (2A;)!! = {2k){2k - 2) ■ ■ ■ (2). 

Note that the result for the scalar relativistic multipoles in d dimensions 

is a remarkably simple generalization of the 3-dimensional result obtained in j^: It is 
enough to replace the explicit 3’s, 5’s etc. appearing in Eq. (B.14b) of [sH hy d, d + 2, etc., 
without changing anything else. In it was also shown that the expansion was in 

3 dimensions the PN expansion of the exact result 

= j dz5f\z) S{j,t +z\y\/c), (2.15a) 

with 6 f\z) = ^ (1 ~ ^ J ^dz 6^°\z) = 1. (2.15b) 

The ratio of Gamma functions appearing in Eq. (j2.15bjl is equal to (2£-|-l)!!/(2^+^£!). Note 
that since the expansion is purely “even” (he., with only even powers of c~^), the time 
argument t -f z|y|/c in (I2.15a|l can be equivalently changed into t — z\y\/c. 

Correspondingly, one can check that the d-dimensional result (im - (imi) is the PN 
expansion of the following simple generalization of the 3-dimensional case: 

= j dz6f\z) S{y,t +z\y\/c), (2.16) 


where we introduced e = d — 3, and 


6f(z) = 


+ - + - 

' 2 ' 2 


) 


r(i)r((-+i + |) 


(1 - 


2\e+i 


dz 6i^\z) 


= 1 . 


(2.17) 


'-1 


Consistently with what happened in Eq. (EH, the kernel 6 f^\z) is simply obtained from 
its 3-dimensional limit by replacing everywhere i hy i + ^ {i.e., 2i hy 2i + d — 3): 

y\z) = Sf^,(z). ( 2 . 18 ) 

Let us mention in passing that the “exact” re-summed expression (ITTBl) can also be directly 
derived from the Fourier-space expression of the d-dimensional Green’s function. 

Finally, having obtained the STF decomposition of the source term in the form 
(ITTTl) . we obtain the corresponding expression of the scalar held (p(x, t). As we pointed out 
above, the remainder term in Eq. (imi) does not contribute to the multipolar expansion 
of the held. Henceforth we shall denote by A4{ip) the multipolar expansion of cp, which is 
therefore given by 

+ CXD / \ £ 

M(v){zc,t) = (Si(«)8ii(x)) , (2.19) 

e=o 

since the terms (□ 5d) give zero when considered outside the compact support of the 

source. In terms of the retarded Green’s function the latter formula becomes 


M(ip)(x,t) 


r ^+00 ^ 

di / dsSL^s) GRet(x, t - s) 

D-n Ld-CXD 


( 2 . 20 ) 


Note that, in view of the retarded nature of the Green function GRet(x, f — s), the integral 
is limited to s < t, and even io s < t — r/c with r = |x|. Equation generalizes what 
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was the basic result for the multipolar expansion of a 3-dimensional inhomogeneous wave 
equation = S, namely 




1 

dvr 


+ O0 

E 


£=0 



St"\t-r/c) 


( 2 . 21 ) 


A common feature of the result and its d-dimensional generalization is that each 

“multipolar wave” of degree i is obtained by an Atuple differentiation, with respect to the 
spatial coordinates, of an elementary spherically symmetric (he., monopolar) retarded solu¬ 
tion; indeed, as mentioned above GRet(x, t — s) depends only on r and t — s. In 3 dimensions 
the elementary spherically symmetric retarded solutions admit a simple expression in terms 
of the multipole moments, namely {t—r/c) jr. By contrast, the d-dimensional analogue 

of each elementary spherically symmetric solution is a more complicated non-local functional 
of Sl{s), which involves an integral over its time argument: ds Sl{s) G^^eti'f'jt — s). 

This non-locality in time in the expression of (f in terms oi Sl comes in addition to the non¬ 
locality in time entering the exact dehnition of the effective source term sf\y,t). 

The former non-locality is evidently related to the fact that the “Huygens principle” holds 
only in d = 3, 5, 7, ■ ■ ■ dimensions. In these special dimensions, the support of the retarded 
Green function Giiet{r,t — s) is concentrated on the past light cone s = t — r/c. On the 
other hand, in other dimensions (and notably in dimensionally-continued complex ones) the 
support of the retarded Green function GRet(^, ^ ~ s) extends over the interior of the past 
light cone: s <t — r/ c. 


III. MULTIPOLE DECOMPOSITION OF THE GRAVITATIONAL FIELD 


A. d-dimensional generalization of the multipolar post-Minkowskian formalism 


27 . lip, applied the general ex- 
UqL n0|, which are themselves 


The calculations of the 3.5PN templates. Refs. j25|, 12: 
pressions of the relativistic multipole moments of Refs 
to be inserted into the (3-dimensional) multipolar post-Minkowskian (MPM) formalism of 
Ref. js^. Let us sketch how one can, in principle, generalize this MPM formalism to arbi¬ 
trary dimensions d. The basic building blocks of the MPM formalism are: 


(i) the parametrization of a general solution of the linearized vacuum Einstein equations in 
harmonic coordinates, say by means of several sequences of irreducible multipole 
moments; 

(ii) the dehnition of an integral operator, called FPDj^g^, which produces, when it is ap¬ 
plied to the non-linear effective MPM source ^ 2 ,..., hn-i) appearing at 

the non-linear iteration, a particular non-linear solution, p/f, of the inhomogeneous 
wave equation Dp/f = 

(iii) the dehnition of a complementary homogeneous solution ql//' {G\ql//' = 0) such that 

h/f = p/Zf -|- ql/Z' satishes the harmonicity condition = 0. 


Given these building blocks, the MPM formalism generates, by iteration, a general so¬ 
lution of the non-linear vacuum Einstein equations as a formal power series, \/—gg^'^ = 










7 ]^^ + Ghi'^ + ■ ■ ■ + + • • •, this solution being parametrized by the arbitrary “seed” 

multipole moments entering the dehnition of the first approximation h^’^. We briefly indi¬ 
cate how the various building blocks can be generalized to arbitrary dimensions d. We have 
in mind here an extension to generic integer dimensions d > 3, before defining a formal 
continuation to complex dimensions. [We consider mainly larger dimensions d > 3 because 
they exhibit generic d-dependent features, while lower integer dimensions, d = 1,2, exhibit 
special phenomena.] 

In the previous section we have discussed the multipole expansion of scalar helds, \3ip = S, 
in arbitrary d. We have seen that the general (retarded) solution outside the source S could 
be parametrized, in any d, by a set of symmetric trace-free (STF) time-dependent tensors 
S'i(t). The situation is somewhat more complicated for other fields, notably the spin-2 field 
relevant for gravity in any d. As we shall discuss in the next subsection, the multipole 
moments needed in a generic d > 3 to parametrize a general gravitational held are more 
complicated than what can be used in d = 3. In d = 3, one can use two independent 
sets of STF tensors, say (the “mass multipole moments”) and Sl (the “spin multipole 
moments”). In a generic d > 3, one has still the analogue of the mass multipole moments, 
he., STF tensors Ml corresponding to a Young tableau made of £ horizontal boxes ( i i i - i i ). 
However, the spin multipole moments must be described by a mixed Young tableau having 
one vertical column of two boxes and ^ — 1 complementary horizontal ones — so that there 
are £ boxes on the upper horizontal row ( |J ^^ I). In addition, one must introduce a 
third type of irreducible representation of the d-dimensional rotation group 0(d), namely a 
mixed Young tableau having two vertical columns of two boxes and ^ — 2 complementary 

^). For instance, when ^ = 2, this new irreducible representation has 

As is well known, this representation 


horizontal ones 


the symmetry of a Weyl tensor in d dimensions; 
does not occur in d < 3. However, all these technical complications will have little impact 
on what we will need to calculate here. Indeed, as discussed below, it will be enough for 
our purpose of unambiguously computing the 3PN-level gravitational radiation emission to 
deal with the simpler mass multipole moments Ml, which admit a uniform treatment in 
any dimension d (actually we shall use a specific definition for what we call the source-type 
mass multipole moments and denote them by II instead of Ml)- 

Let us turn to the generalization of the integral operator FP In d = 3, the precise 

dehnition of this operator was the following. Consider a typical non-linear source generated 
by the MPM iteration, e.g. N 2 = N 2 {hi, hi) ~ (dhi)"^ -|- hid'^hi, in which hi is represented 
by its multipole expansion. One formally assumes that the multipole expansion of hi ~ 
'Y^d[ML{t — r/c)/r] contains a finite number of multipoles. This ensures that the non¬ 
linear source N 2 {hi) is a hnite sum of terms of the form hKF{t — r/c)/r'?, with angular 
factor hk = STF (n'l ■ ■ ■n*''), n* = jr. We can further expand Fit — rjc) in powers of 
r/c and get ^ 2 (^ 1 ) as a sum of terms ~ hxF{t)/rP. Though this multipole expansion of 
N 2 {hi) is only physically relevant in the region outside the source, say r > a, in the MPM 
formalism we always mathematically extend its dehnition (by real, analytic continuation in 
r) down to r = 0. Then this formal construction, h-MPM = G hi + h 2 + ■■■, valid by real 
analytic continuation for any r > 0, is identified with the multipolar expansion, say Ai{h), 
of the physical held h. While the physical h takes diherent expressions inside (r < a) and 
outside (r > a) the source, the object M.{h) = hupM is mathematically dehned everywhere 
(except at r = 0) by the same formal expression, but is physically correct only when r > a 
(see for the notation and further discussion). 
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To deal with the singular behavior near r = 0 of the non-linear MPM source terms, e.g. 
N 2 {hi) ~ nKF{t)/r^, one introduces a complex number B and considers the action of the 
retarded Green operator onto the product of the source by a “regularization” factor (r/ro)'®, 
say 


p{d—3) 
^2 


(s) ^ ni 



(3.1) 


The length scale tq represents an arbitrary dimensionful parameter serving the purpose of 
adimensionalizing the above regularization factor. It was shown in Ref. 5^ that the integral 
F 2 {B), Eq. (jSHI), is convergent when the real part of B is large enough, and that F 2 {B), 
considered as a function of the complex number B, is a meromorphic function of R, which 
has in general (simple) poles aX B = O,'^ coming from the singular behavior of the integrand 
N 2 {hi) near r = 0. [One formally assumes that the multipole moments are time-independent 
before some instant —T, and at the end of the calculation the limit T —>■ -|-cx) is taken.] 
Therefore, the Laurent expansion of F 2 {B), near R = 0, is of the form 


^^(d-3)(^) ^ ^ ^ Ci{^,t)B + O (R2) . (3.2) 

One then dehnes, when d = 3, the finite part at B = 0 of n^l^^N 2 {hi), denoted 
FP □j^g^.A^ 2 (hi), as the term Go(x, t) in the Laurent expansion of F 2 {B). One proves that 
Go(x, t) satishes the equation nOo = A^ 2 (^i) and uses it as the “particular” second-order 
contribution pi^'^ to the second-order metric . Let us not spend time on the construc¬ 
tion of the additional homogeneous contribution qtfi' necessary to satisfy the harmonic- 
ity condition d„ {pifi' -|- g^^) = 0 [an example of construction of such contribution will be 
given in ()3.41j) below]. Having so constructed (in d = 3) the second-order term in the 
MPM expansion of the external metric, = p'fi' + , one continues the iteration by 

considering the next order inhomogeneous equation nh^ = N^{hi,h 2 ) and introducing 
FfiB) = \fir/ro)^Nfihi^h 2 )\- The singular behavior near r = 0 of is more com¬ 
plicated (it contains logarithms of r), and, as a consequence, one finds that FfiB) is still 
meromorphic in the complex B plane, but will contain double poles at R = 0. Again, one 
dehnes p^ = FP D^g^. A 3 as the coefficient of the zero-th power of B in the Laurent expansion 
of FfiB) when R ^ 0. 

Having recalled the dehnition and properties of the operation FP in the 3- 
dimensional MPM formalism, let us sketch what changes when working in d dimensions. 
Let us start with the “seed” linearized metric hfi'. As we see in Eq. fl2.19jl . and will see 
below with more details for the tensorial analogue of the scalar multipole expansion, the 
multipole expansion hi is of the form hi ~ ^ [Mi,(f)h(x)]. Though one cannot write, 

in arbitrary d, a simple, closed-form expression for the object D^g^. [Mi(t)d(x)], it is enough 
to write down its expansion when r —> 0 (which is in fact the same as its PN expansion). 
Modulo regular terms near the origin, this expansion is obtained as 


□ 


-1 

Ret 


MLit)6{x.) 


-|- regular terms. 


ML(t)S{x) 


(3.3) 


^ Actually, it was shown in mm that F 2 {B) happens to have no pole when B ^ 0, due to the particular 
structure of the quadratic-order interaction. 
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Using A“^(5(x) oc A“^5(x) oc etc. we see that the 3-dimensional form of the 

expansion of hi near r = 0 (after taking into account the expansion of the retardation r/c), 
takes in d dimensions the form 


hi ~ ^ 


nKF{t) 

rfP+S ’ 


(3.4) 


where uk = STF[nj^ ■ ■ p is a (relative) integer, and e = d —3. Inserting this expansion 
in the second-order source N 2 {hi) ~ dhidhi + hi&^hi yields 


N 2 (hi) ~ 


nKpit) 

' rpp-\-‘2s 


(3.5) 


At this stage, one could consider without inserting a factor (r/rg)'®, by using the an¬ 

alytic continuation in d. However, to ensure continuity with what was done in 3 dimensions, 
it is better to insert this factor and to consider 


= □ 


-1 

Ret 



B 

Ni (fti) 


(3.6) 


The main difference between (Eini) and its 3-dimensional analogue (EH) concerns the mero- 
morphic structure of F 2 {B). Indeed, in view of the shift by +2e of the integer exponent p 
in (|3.5p . and of the presence of in the d-dimensional volume element d'^x = r‘^^'^dr dVL 2 +£, 
one easily sees that the (simple) poles in F 2 {B) that were located at H = 0 when d = 3 are 
no longer located at H = 0 when d 7 ^ 3, but are shifted ai B = 2e — e = e. Alternatively, 
this can be explicitly verihed by using the expansion = A“^ -|- -(-■■■ (plus a 

regular kernel), and the formula A“^r“ = /[{a -f- 2 ){a + d)] where the pole at a = —d 

is the only one which comes from the ultraviolet (UV) behavior r —>■ 0. As a consequence, 
the expansion dH is now modihed to 

f(^\b) = + C'J‘'^(x,t) + c[^\^,t)B + O [B^] . (3.7) 

B — £ 

This expansion, and its analogues considered below, is considered for e and B both small (so 
that the expansion in powers of B makes sense), but without assuming any relative ordering 
between the smallness of B and that of e. One should neither re-expand {B — £)~^ in powers 
oi Bje nor in powers of e/B. 

Having in hand the above structure, one then defines the d-dimensional generalization of 
the hnite part of N 2 {hi) as the coefficient of B^ in Eq. ()3.7D . namely Cq^\ We denote such 
a finite part by 


or, more simply, by 


FP □Re\ 


r^N2 {hi] 
where r = 


FPnoiAfRftO = FPK 


Id) 


C'(S‘^^(x, t ), 

(3.8a) 

r 

(3.8b) 

ro 

III 

(3.9) 


Note the subtlety that the expansion (Eli) is neither a Laurent expansion in powers of B — e, 
nor a Laurent expansion in powers of B. After subtracting the shifted pole terms cx {B—£)~^, 
one expands the remainder in a regular Taylor series in powers of B. The interest of this 
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specific definition is the fact that it ensnres that (x, t) is an exact solntion of the eqnation 
we initially wanted to solve, namely 


n^d)c(d) ^ N2{hi). (3.10) 

Indeed, by its mere definition dsini), one has = {r/ tq)^N 2 {hi). Comparing this 

result (which has no pole) to the application of □ to ()3.7|1 . we first see that the pole part 
must be a homogeneous solution, = 0. Then, identifying the successive powers of B 

(using = 1 + Blnr + ■ ■ ■), yields = N 2 , = ln(r/ro)A^ 2 , and so on. 

Another useful property of the d-modified definition (ESI) is that it automatically ensures 
the continuity between d —>■ 3 and d = 3. Indeed, the shift in the location of the pole in 
(JUl) was made to “follow” the pole that existed at i? = 0 when d = 3. Therefore we have 
lim^^s = C-i, and similarly lim^^aCo'^^ = Cq, etc., where the RHSs are those defined 
in Eq. dSSD when d = 3. 

The extension of the iteration to higher non-linear orders introduces a new subtlety. 
Indeed, let us look more precisely at the structure of the second-order contribution to the 
metric, h 2 = P 2 + (I 2 where, as we said, the particular solution p 2 is dehned by the modihed 
FP process: p 2 = FP □^g^A" 2 (hi), and where q 2 is a complementary homogeneous solution. 
Most of the terms in the integrand N 2 introduce no poles, and, for them, we simply hud 
a structure of the type ~ simplicity, we henceforth suppress angular 

factors). Let us now consider the terms in N 2 that generate poles oc {B — e)~^ in F 2 {B).^ 
We know that such terms introduce, when d = 3, some logarithms of the radial variable 
r. When d ^ 3, they no longer introduce logarithms but they introduce a further technical 
complication. Indeed, let us look at a typical example, namely a dangerous term in F 2 {B) 
of the form F 2 ^°^^\b) = /\-i( 7 >^- 3 - 2 £^ Suppressing for simplicity a factor (B — 1 — 2e)~^ 
which is jointly analytic in B and e near B = 0 and £ = 0 respectively and therefore creates 
no problem, we have essentially F 2 ^°^^\b) = [B — e)~^r^~^~^^. According to Eq. (13.7j) the 
pole part of that we must subtract is, for instance, obtained by multiplying hjB — e 

and then taking the limit B ^ e (and not i? —0). This pole part is therefore given by 
[B — e)~^r~^~^. The hnite part of F^°''^\b) is then obtained by subtracting the pole part 
and taking the limit R ^ 0 ; this yields 

p(pole) _ pp p,^{pole) ^ ^ 

The subtlety is that there seems to appear poles in However, the residue of the pole 
vanishes, since the limit £ —0 of dTTTl) is finite and generates the logarithm that we know 
to exist in d = 3, ^ 2 *^°^*^^ Inr/r. If we do not take the limit e —>■ 0, we must keep the 
structure (IXTTD and see what it generates at the next, cubic, order of iteration. In addition, 
we must also add the complementary solution q 2 needed to satisfy the harmonicity condition, 
= 0- calculation of q 2 could be done in d = 3 without encountering 

poles (see Ref. jl^), it will clearly not create problems in d 7 ^ 3 apart from the fact that 


® As said above, these terms actually cancel among themselves because of the particular structure of A^ 2 - 
However, similar terms appear at higher iteration orders, and their general structure is simpler to describe 
if we start our induction reasoning at the quadratically non-linear level. 


12 





q 2 , being a homogeneous solution \I\q 2 = 0, will behave near r = 0 essentially like hi, he., 
<?2 ~ (which differs from most of the terms of p 2 which were 

Summarizing so far, the second-order MPM iteration ^2 has a structure, near r = 0, of 
the symbolic form 

h 2 ~ Cl (e) r-P-^ + C 2 (e) ^ , (3.12) 

where the Cj(e:)’s are analytic at e = 0, and where we explicitly separated the semi-singular 
structure in e. When inserting the structure dSH into N^^hi, ^ 2 ), one finds that the singular 
behavior of N 3 near r = 0 can generate several types of singularities in B and £. There are 
simple poles oc {B — e)~^ and simple poles oc {B — 2e)~^, which are natural generalizations of 
the structures that generated simple poles in F 2 {B). When looking at the effect of the more 
complicated structure given by the third term on the RHS of (Em, one finds that it is best 
described as generating some “quasi-double poles”, namely terms oc {B — e)~^{B — 2e)~^. 
The point is that if one were to expand this term in simple poles with respect to B, namely 


{B-e){B-2e) e{B-2e) e{B-e)' ^ 

it would seem to involve poles in 1/e. However, all such poles are “spurious” because the 
source of the trouble which is the last term in ()3.12|1 had a finite limit as £ —>■ 0, and because 
one can easily see that, in our above-defined MPM algorithm, source terms having a finite 
limit as e —0 generate solutions having also a finite limit as e —0. 

Finally we find, by induction, that at each iteration order n one has the structure 

+ ■ ■ ■ + dn{e) , (3.14) 

where the coefficients di{e) might individually have (simple or multiple) poles in e, e.g. 
di{e) = Ci{e)/e^, but which always compensate each other in the complete sum Then we 
obtain that the integral 




r®iV„ (/ii, ■ ■ ■ 


(3.15) 


will have an expansion, near i? = 0, of the generic form® 

= E F) ■ (3-16) 

The “quasi-multiple poles” which constitute the first term on the RHS have k <n — 1 and 
1 E Q'i < — 1 • As we have seen in flTT^ and dsnii the poles in 1/e are in fact spurious, as 

they have a residue which is always zero. [We assume here that the seed multipole moments 
are regular as e —>■ 0.] So, when writing the result in the form of duni), we note that the 
coefficients C'q'^\ etc., are all regular when d —>■ 3. One then defines the d-dimensional 


6 


Our notation is a little bit over simplified, since the coefficients c/// depend in fact on a set of integers 
{qi, ■ ■ ■ , qk\. Also we do not indicate the obvious dependence of the coefficients on the iteration order n. 
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generalization of the finite part of as being the coefficient of in the expansion (ITTBD : 

= (3.17) 

This coefficient is regnlar when e: = cf — 3 —>■ 0, thongh it contains apparently singnlar terms 
of the type of the last term on the RHS of ()3.12|1 . Moreover, nsing the same reasoning as 
above, one finds that it satisfies the needed resnlt: = Nn- Note finally that, when 

e —0, the qnasi-mnltiple poles in (j3.1fi|l merge together to form the mnltiple poles oc B~^, 
with fc < n — 1, that were fonnd to exist in d = 3 jl^. On the other hand, when e 7 ^ 0, 
the poles form a “line” of simple poles located ed, B = e, B = 2e, ■ ■ ■, B = {n — l)e. 
However, it is better not to decompose the product of simple poles entering (I3.1b|l in sum of 
separate simple poles, because this decomposition would, as in Eq. dm, introduce spurious 
singularities (x e~K 

The main practical outcome of the present subsection is the modihed definition of the 
operation FP when working in d 7 ^ 3: namely as the coefficient of B^ in an expansion 
of the type (IXTBD where, after separating the shifted poles at B — e, ■■■, B — {n — 1)£, one 
expands the remainder in a Taylor series in powers of B. Note that a simple consequence 
of this dehnition is that, for instance, a term of the form B/{B — qe) in Dj^gt gives 

rise to a finite part equal to 1. Indeed, 


FP 


B 


B — qe 


= FP 


B — qe + qe 
B — qe 


= FP 

\J" +1I 


B — qe 


= 1 


(3.18) 


One might have been afraid that a term of this type, B/{B — qe), could have been ambiguous 
because, ultimately, we are sending both B and e towards zero without fixing an ordering 
between the two limits. However, in the present d-dimensional generalization of the MPM 
formalism, everything is precisely defined and unambiguous. 


B. Multipolar decomposition of the gravitational field in d dimensions 


We now sketch the d-dimensional generalization of the results concerning the matching 
between the MPM exterior metric and the inner field of a general post-Newtonian (PN) 
matter system. To start with, we consider the case of a smooth matter distribution, and 
will later allow the matter stress-energy tensor to tend to a distribution localized on some 
world-lines. The next subsection will be devoted to the d-dimensional definition of the 
source multipole moments. The investigations of this and the next subsection are based 
on the works Id^ 1 ^ which derived the expressions of the source multipole moments 
of a general PN source, up to any PN order (in 3 dimensions). Early derivations of the 
relativistic moments, valid up to IPN order, can be found in Refs, j^, . 

We look for a solution, in the form of a PN expansion, of the d-dimensional Einstein field 
equations. As before we choose some harmonic coordinates, which means that = 0 

where the so-called “gothic” metric deviation reads = \/—g , where g is the 

determinant and g^'" the inverse of the usual covariant metric gpa- Then the Einstein field 
equations, relaxed by the harmonic coordinate conditions, take the form of some “scalar” 
wave equations, similar to (EH), for each of the components of , 




IGttG 


■ T 


fll/ 


(3.19) 
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where □ denotes the d-dimensional flat space-time wave operator, and G is the d-dimensional 
Newton constant related to the usual Newton constant Gat in 3 dimensions by Eq. (113) 
below. The main contribution we shall add in the present subsection, with respect to our 
investigation of the scalar wave equation in Sec. im is how to deal with the crucial non-linear 
gravitational source term in the Einstein field equations, which makes the RHS of Eq. (ixm 
to have a support which is spatially non-compact. The RHS of involves what can be 

called the total stress-energy pseudo tensor of the non-gravitational and gravitational fields, 
given by 


= + (3.20) 

where is the matter stress-energy tensor, and the second term represents the gravita¬ 
tional stress-energy distribution, which can be expanded into non-linearities according to^ 

AM- = Af (h, h) + Ar(h, h, h) + • • ■ , (3.21) 

where the quadratic, cubic, etc., pieces admit symbolic structures such as A 2 rsj hd‘^h-\-dh dh 
and A 3 ~ h 9/1 dh. 

The solution of the held equations we consider in this subsection will be smooth 
and valid everywhere, inside as well as outside the matter source localized in the domain 
r < a. Inside the source, or more generally inside the source’s near zone (r -C A, where A is 
the wavelength of the emitted radiation), will admit a PN expansion, denoted here as 

. On the other hand, in the exterior of the source, r > a, will admit a multipolar 
expansion, solution of the vacuum field equations outside the source, and decomposed into 
(d-dimensional) irreducible spherical harmonics. As usual, the definition of the multipole 
expansion is extended by real analytic continuation in r to any value r > 0. It will be 
necessary to introduce the special notation Ad(hM-) |;o mean the multipole expansion of 
As we already mentioned, the multipole expansion in the present formalism is given 
by the MPM metric of Sec. nnn which is therefore in the form of a formal infinite post- 
Minkowskian series up to any order n, 

M.{h^ ) = ^MPM • (3.22) 

As mentioned above, though the identification ()3.22j) is only physically meaningful in the 
exterior domain r > a, it can be mathematically extended down to any r > 0 by real analytic 
continuation in r. 

In this subsection we shall show how to relate in d dimensions the multipolar expansion 
(Km to the properties of the matter source, in the case of a PN source (he., one which 
is located deep inside its own near zone, a -C A). Actually, the derivation below will be a 
simple d-dimensional adaptation of the proof given in the case of 3 dimensions in Ref. |4^ 
(see notably the Appendix A there). 

The heart of the method is to show that one can deal with the presence of non-compact- 
support source terms on the RHS of the held equation by considering a certain 

quantity which satisfies a wave equation whose source does have a compact support. 


^ In the MPM formalism of Sec. IIII Al we used N 2 {hi) = A 2 (/ii,/ii), N^(hi,h 2 ) = -f 

A 2 (/ii, / 12 ) -t A 2 (/i 2 , hi), and so on. 
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and thus, whose multipolar expansion can be computed by using the results of Sec. HT] (for 
each space-time component fiu). This quantity is defined by 


^ _ pp □-! 




(3.23) 


The second term in (j3.23|l . that we thus subtract from in order to dehne this quantity, 
involves the hnite part operation FP in d-dimensions which has been dehned in the previous 
subsection (uraj. It contains the regularization factor = (r/rg)^. The use of the 
operator FP is consistent with Sec. 1111 Al because it acts on the multipole expansion of 
the non-linear source term which is in fact identical to the formal post-Minkowskian 

(MPM) inhnite series cf. Eq. (I3.22|l . The meaning of the last term on the RHS of 

(non is that FPDpgj. is to be applied to each term of the MPM expansion of r^AA{A^‘'), 
and that we then consider the formal summation of this MPM series. 

Equation appears to be the difference® between the solution of the field equation 

and the contribution coming only from the non-linear terms in the exterior of the 
compact-support source (and then analytically continued down to r = 0). Since is the 
retarded integral of the pseudo tensor and since the multipole expansion of the matter 
tensor is formally zero: = 0 (because has a compact support), we can rewrite 

fj3.23|l as 


= 


Ibvr G 


r-i-l _ pp 

'-'Ret 




(3.24) 


Next, we remark that the hrst term in is regular within the source (for r < a), and 

that we can therefore add to it the same FP procedure as in the second term, without 
changing its value — because for regular sources, the operator FP simply gives back 
the usual retarded integral. Thus, 




IGttG 1 

-^ FP Dhet 

c4 B ^ * 




A4{r 




(3.25) 


As we said, the multipole moment formalism we are using is defined for general smooth 
matter distributions [say G C'°°(M^)] with compact support. Hence, is regular inside 
the source, and is a perfectly well-dehned object. Only when general formulas for 

the multipole moments are in hand, shall we apply them to point particles (in Sec. and 
then shall we need a self-field regularization scheme to cure the divergencies induced by the 
point-particle model. Of course the FP procedure used here should be carefully distinguished 
from the self-held regularization. 

The point is that in the form given by Eq. (nnHU . appears now as the retarded 
integral of a source with compact support (limited to r < a). This follows from the fact 
that agrees numerically with its own multipole expansion in the exterior of the 

source, for r > a. Hence we are allowed to use the end results of Sec. H which applied to 
compact-support sources (and those results can evidently be “uniformly” applied to all the 
components of A^*^). From and we obtain 


M(A'^") 


IGttG 


+ 00 

E 


i=0 



b+oo 




s) G 


Ret 


.X, t 



(3.26) 


® This “difference” has of course nothing to do with the difference between the dimensional and Hadamard 
regularizations that we consider in Sec. El 
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(3.27) 


where the multipole-moment functions read 


= FP 


■iVlyl"?! / dz5f(z) (y,i + z|y|/c). 




A difference with the multipole moments considered in Sec. O is the presence of the FP 
process with regularization factor lyl-® = |y/'ro|'®. We shall see later the crucial role played 
by this FP process. 

Now we deal with the integrand appearing in the multipole moments (I3.27j) . 

following the same argument as in Ref. . Such an integrand has a compact support limited 
to r < a, so we see that in the case of a PN source, for which a A, we can replace it with its 
formal PN expansion, because precisely the PN source is conhned within the source’s near 
zone. Hence the PN-expanded moments will be generated by the PN-expanded integrand 
where we denote the formal PN expansion with an overline. 

Let us now show that the second term, M. , gives no contribution to the PN moments. 
We know that the structure of this term (in d = 3 -f £ dimensions) reads 


where p and q are relative integers. This follows from Eq. ()3.14|1 above. The argument 
showing the vanishing of the term involving is that any term of the type ()3.28|1 

in the moment will ultimately give [after taking the PN expansion like in ()2.14|l ] a spatial 
integral of the type J d‘^:s.nK ^ say (times some function of time), which we know to 

be exactly zero by analytic continuation in B. Therefore, following this argument, which is in 
fact the same in d dimensions as in 3 dimensions, we are led in fine to a PN multipole moment 
which is simply generated by the PN expansion of the (non-compact-support) pseudo tensor, 
T^'^. Hence, we write our result as 


= j d'^y\y\^yLr^[]{y,t), (3.29) 

where the ^-dependent integrand takes the form of the following PN expansion, 

^[i]{y,t) = j dz6f\z)T>^''{y,t + z\y\/c) 

+00 y , I p. \ 2k 

k=0 ^ ' 

The PN coefficients have been given in ()2.9bfl . Note that the final result in (nom 
combines two separate PN expansions: (i) a PN expansion of the type ()2.13j) (already 
indicated by an over line notation), and (ii) the usual PN expansion of The presence of 
these PN expansions is crucial to the meaning and validity of the final expression in ()3.3()j) . 
Finally, note that our use (in the proof above) of the vanishing of the spatial integrals 
J d’^'x.nK implies that we have transformed the role of the factor |y|^ from that 

of regularizing integrals that are singular at r = 0, into that of regularizing integrals that 
are singular at r = oo. Thereby, in the final result (USED, the FP procedure is used as a 
regularization of the boundary at infinity of the integral, which would otherwise be divergent 
because of the multipolar factor ~ |y|^ multiplying the non-compact-support (and PN- 
expanded) 
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C. Symmetric-trace-free source multipole moments in d dimensions 


In Eq. we have represented the quantity A^'^, Eq. (Em, in the form of an inhnite 

superposition of scalar multipolar waves, say where we associate to any function of 

time {t) a corresponding spherically symmetric retarded wave denoted® 


^+oo 






Ret 


'x,« - s), 


(3.31) 


in which the tensor indices play the role of simple “spectators”. This expansion is not 
yet a genuine irreducible tensorial multipole expansion. To transform Eq. in a ten¬ 

sor multipole expansion, we need to decompose each “elementary wave” in irreducible 
representations of the d-dimensional rotation group 0{d). As each (undifferentiated) ele¬ 
mentary wave Eq. (Id.dlj) . is spherically symmetric, the problem of decomposing 

iF'^’^{r,t) in irreducible components is reduced to the purely algebraic problem of decom¬ 
posing its “source” whose expression is given by Eqs. (Id.29j) - (ld.d()|l . in irreducible 

representations of 0 {d). 

Let us consider in turn the various components of The time-time component 

is already put in irreducible form because it is STF with respect to the multi-index 
L. In the language of Young tableaux [for 0(d)], the STF-£ representation carried by 
is denoted by ^ horizontal boxes i i i i i . The time-space component JF®* is, algebraically, 
the product of an irreducible vector representation Y* and of an irreducible STF-£ one Tj,. 
In Young tableaux terms, this corresponds to the product nx i i i - n . In any dimension d, 
this product gives rise to three irreducible representations: the STF-(£ + 1) one i i i i -1 i , 
the STF-(f' — 1) one i i - i i , and a mixed Young tableau representation |J ^^ L The hrst 
two representations are easily understood as corresponding to the STF projection of 

the product I^T^, and the contraction K^aL-i- It is more intricate to write explicitly the 
mixed-Young-tab lean representation contained in ViTi. When d = 3, it was convenient to 
use the Levi-Civita antisymmetric tensor Sijk to dualize the two antisymmetric indices in the 
mixed Young tableau (he., the two vertical boxes) and replace them by a single vector index. 
However, when considering a generic value for the dimension d (initially taken as an integer, 
and then formally continued to arbitrary complex values), one is not allowed to use the 
specihcally 3-dimensional tensor Sijk- When further considering the space-space component 
one is facing the algebraic problem of decomposing the product (nxni..n,x [ i i - i i . Again 
the irreducible decomposition of this product contains both relatively simple symmetric 
representations, such as i i i i i n , and more involved mixed-Young-tableau ones of the 
IJ II I or “Weyl” —type.^® It is quite possible to write this decomposition 


■ spin 


in any (integer) dimension d, by using mixed-Young-tableau projectors instead of the Sijk 
duality operations used in d = 3. However, for the present work we can simplify our task 


® In 3 dimensions we recover t) = {t — rlc)lr. Recall that in any dimension the Green function 

GR,et(x, t) is in fact a function of r = |x| and t. 

Note in passing that in the irreducible decompositions of AjJ* and in 3 dimensions, the Levi-Civita 
tensors always appear in pairs, and that the products ee which appear can always be entirely expressed 
in terms of Kronecker deltas: oc = SaS^Sc^. 
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and simply ignore all the mixed tableanx that appear in the irredncible decomposition of 
and 

Indeed, becanse of their irreducible character the symmetric Young tableaux rrr 


i: 


never mingle with the mixed-symmetry ones when doing linear operations, as is done when 
working with a linearized solution of Einstein’s equations. (This is easily checked when 
looking, for instance, at the derivation of the multipolar expansion of linearized gravity 
given in jl^.) Moreover, when considering the gauge-invariant content of (he., modulo 
a linearized gauge transformation), the symmetric representations will hnally give rise to the 
so-called mass-type multipole moments, say II, while the mixed one will correspond to what 
are called (in d = 3) the spin-type or current-type multipole moments, say Jl. The aim of 
the present work is to cure the ambiguities, linked to logarithmic divergencies of integrals 
calculated in Hadamard’s regularization, which appeared at 3PN order in the calculations 
of Refs. 


However, the works j26|, |23,12^ found that the only multipole moment 
which introduce ambiguities (in d = 3) at the 3PN order is the mass quadrupole moment lij. 
Therefore, for our purpose, it is enough to derive general d-dependent formulas for the mass 
multipole moments II (besides the quadrupole lij we shall also need to consider below the 
mass dipole moment Jj). We do not need to consider the dehnition of the current moments 
Jl outside of d = 3. In view of what we said above, it is therefore enough to consider only 
the symmetric ( i i i -1 i ) pieces in the irreducible decomposition of dq?-’ and We then 
write 


•rOO 

L 

riO 

L 


— Rl 1 


-iL 

+ STF STF 


h:’i) -|- “mixed tableaux” 


+ 


L i] 

“mixed tableaux” 




r(“2) 


jL-l + ^iH^jH-iRL-2 


+ S^jVl 


(3.32a) 

(3.32b) 

(3.32c) 


where the angular brackets surrounding indices refer to the STF projection, and where the 
tensors Rl, dl+i, ^ ^TF in their indices (recall our notation 

for multi-indices: L = ii ■■■ ie, L + 1 = ii ■ ■ -ii+i, etc.). Furthermore, we shall need below 
the inverse of Eqs. (im^ . he., the expressions of these tensors in terms of the jF’s. These 
are 


Rl 

'TOO 
— L 

(3.33a) 

^L+l 

~ ^L) ’ 

(3.33b) 

rpi — ) 

(.(2(- + d-4) 

(£-hd-3)(2£ + d-2) 

(3.33c) 


^(*^+2*£+l 

(3.33d) 

Y* 

2rf<(2(’ + rf-4) „ 

(d - 2)(£ + d - 2)(2£ + d) L ’ 

(3.33e) 

^(-2) 
^ L-2 

£(£-l)(2 £ + d-6 ) 

(£ + d - 3)(£ + d - 4)(2£ + d - 2) “^^-2 ’ 

(3.33f) 

Vl 

^ ^aa 

- • 

(3.33g) 


The next step towards the dehnition of the STF source moments is to take into account 
the effect of the harmonicity conditions = 0) on the multipolar expansion (I3.2bj) . 


19 














which we henceforth write with the help of the shorthand notation (ICTl) as 




AG 


+00 

£=0 




dT.:F- 


■flU 




(3.34) 


The latter tensor is not divergence free in the full non-linear theory. Indeed, by 

using the same method as the one employed in 3 dimensions and which resulted in Eqs. (4.5)- 
(4.6) of 0, i.e., by using the explicit expressions of the multipole moments, 

we can derive the following relation 






2 i + d 


^aL = 0^ , 


(3.35) 


where the dots mean the time differentiation, and where the new “multipole-moment” func¬ 
tion is given by 

G^{t) = YV j d‘^yB\yfyL^T^;^iy,t). (3.36) 

Recall that r|),“(y,f) is dehned by Eq. (I3.3()j) above.Notice that Eq. ()3.36|1 . involves 

an explicit factor B in its integrand. Hence, by the properties of the analytic continuation 
in B, Qf{t) depends in fact only on the behavior at the boundary of the integral at inhnity, 
|y| ^ -|-cx) (hence, for instance, there is no analog of Gf{t) in linearized gravity). The 
consequence of Eqs. (JSSSD-dlSED is that the harmonicity conditions can be expressed by 
saying that the divergence of the multipole expansion A^(A^^) reads 


+ 00 




c’ -- i\ 
£=0 

/‘+00 


dxGi , 


where Q^ir^t) = -dvr / dsQf^{s) - s) 


(3.37a) 

(3.37b) 


Next we decompose the components of the function Gf in STF guise, which means 

Gl = Pl, (3.38a) 

Gl = Qit'’ + + “mixed tableaux” , (3.38b) 


The result is easily checked once we remember from Eq. that the kernel function 

in d dimensions is equal to 6 ^^e{z). So, the same method as in 3 dimensions applies with the 
simple replacement £ —>■ i + ^. In particular one uses in this derivation 



dz'^ 



-{2t + e + 3)zS^/\z ), 


-(27 + e + 3)(2£ + £ + l) S^/\z) - S, 
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where the tensors Pl, Q^l^i and are all STF, together with the 

inverse formnlas 

Pl = a, 

(3.39a) 

/n(+) _ Pde+i 

Vi+i — yv) ’ 

(3.39b) 

l{2t+d-4) - 

{i + d - 3){2i + d - 2y'^^-^' 

(3.39c) 


The relations (HSZD can then be re-stated as the following constraint eqnations linking the 
STF tensors (to simplify, we set c = 1 for a while), 


Rl = + + 


d + i-2 




= 


t(-) = 


{^ + l){d + 2^-2) 

p(+) I (n {d-2){d + i-2){d + 2i) --qo) 

^ 2id{d + 2i-4){d + 2i-2)^^ + 
{d-2){i+l) ^(0) ^ + 


(3.40a) 


Qh + 


Orl 




Onr aim is now to obtain the linearized mnltipolar solntion, which is at once solntion of 
the sonrce-free eqnations and divergenceless, and which will be exactly eqnal to the linearized 
metric of the MPM formalism. To this end, we introdnce the object given by the 
following mnltipole expansion in d dimensions [recall onr notation (j3.31 jl ]. 






11 

Qi 


-4G 

-AG 

-AG 


-fP + a ( jp + JJqM - 


2d 


(-1-) , 3d -l- 1 1 


-IQ^ + 


2d 


-E 


l>2 


l\ 


dl-lPiL-l 


l>2 


(-) 


i(-) 


^ £1 {^^ij^L-lQi-i Q dL-2{iQj)i_2 


+dL-2 ( PijL-2 + ^ Q‘^L-2 ~ 


7£ + 3d - 6 

ij+v){2rTJ^ 


~(-) 

Q ijL-2 


(3.41a) 

(3.41b) 


(3.41c) 


Here the integral signs refer to a time anti-derivative, e.g. f P{r,t) = j^_^dTP{r,T), 

If Qi'^P) = /!oo/Too ^0- The object , which is given here modnlo the mixed 

tableanx corresponding notably to the spin-type contribntions, exactly corresponds to the 
so-called “harmonicity algorithm” of Ref. (in the slightly modihed version of it proposed 
in Eqs. (2.12) of [s^; notice that the latter eqnations are valid in any d). The properties 
of are that = 0 and [A1(A^^) + qi^] = 0, as one can easily verify by direct 

calcnlation.^^ One now introdnces the object 


hf = M(A^^) + gr- 


(3.42) 


Remember the presence of time anti-derivatives in (ITITll . In the present formalism the metric is past¬ 
stationary, and from this one can show that the functions involved which need to be time integrated are in 
fact zero in the past, before the instant —T, so that there is no problem in defining these anti-derivatives. 
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As in one easily checks that dehnes a linearized multipolar metric (in harmonic 
coordinates), which generates, by MPM iteration, the full metric The source 

multipole moments are then dehned as those which parametrize h^’^. The “main” multipole 
moments will be those which parametrize a specihc piece of the linearized metric sometimes 
referred to as the “canonical” metric and which was introduced long ago in Ref. jH^. The 
canonical metric, say ^, is separately divergenceless, and differs from by a linearized 
gauge transformation, with gauge vector say - 01 , 


hT = Cm + • 

It explicitly reads (still consistently omitting the mixed tableaux) 


where the tilde over the objects II means, as in (IH3T1) . 


^+oo 


hir, t) = -dvr / ds Il{s) GRet(x, t - s). 


(3.43) 


+ 0O / 

£=0 

+0O , Y 

(3.44a) 

£=1 

+00 / Y 

(3.44b) 

4 G ^ ^L-2 AjL-2 , 

e.=2 

(3.44c) 


(3.45) 


Such expressions clearly yield a precise dehnition of the mass-type STF multipole moments 
I Lit) in d dimensions (and the mixed tableaux could be used to dehne some other, “spin- 
type” and “Weyl-type”, moments). We need now to relate the moments II entering ()3.44ll 
to the STF tensors which were used in the STF decomposition ()3.32j) of the function 

To this end it is most convenient to consider the gauge invariant linearized curvature (in 
d dimensions) associated with the metric deviation , in order to eliminate the irrelevant 
linearized gauge transformation in Eq. The component OiOj of the curvature, in 

terms of the gothic metric deviation, reads 

-h 2dod(ih{^^ + dlhd^ . (3.46) 

Since and differ by a gauge transformation, Eq. (j3.43|l . we necessarily have 

^ 0 i 0 j[di ~ ^cani] = 0, which immediately gives us (looking at the particular term pro¬ 
portional to the double gradient didj) the looked-for expression of the moment Ii as being 


h = 


Rj 


d 


:V, 


2 (c^-l) ^(-) 


d-2 ^ c(£ + l)(d-2) 


T)-> + 


d — 1 


c^{^ + l){^ + 2){d-2) 


U 


(-2) 


2(^-3) ^(-) 

{^ + l){d-2y^ ’ 


(3.47) 

where the explicit powers of c have now been restored. We hnd that II is given in terms of 
the STF tensors parametrizing our original mult ip ole-moment function {t) and dehned 
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by Eqs. (ESS), and also, in the last term of Eq. (ITiTli . of the “harmonicity” fnnction 
Giit) given by (ld.df)|l . Note that the last term of Eq. ()d.47j] involves a factor {d — 3) and 
therefore is absent in the 3-dimensional formalism of Since involves also a factor B 
in its integrand, we see that the contribntion indnced in the moments by this term will be 
proportional to B{d — 3)] we shall see that snch a contribntion is actnally zero. 

Once we have obtained the moment Iiit), it is better to express it back in terms of the 
original fnnction since we know its relation to the psendo tensor of the sonrce, given 

by Eqs. Using the inverse relations and dssni), we get 



c%d + i-2){d + 2i)^^^ 


4{d + 2i-2) 


c\d + £ - l)(d + £ - 2)(d + 2£ + 2) 


2 {d + 2l-2) 




(3.48) 


If we fnrther introdnce the following “sonrce-rooted” qnantities 


2 (d - 2)r“ + U* 


(3.49a) 



(3.49b) 

(3.49c) 


we can then write the moment into the following more explicit form 



c2(d + £ - 2)(d + 2£) ^d^+i](y 


4(d + 2£-2) 


^iLSi[£+i](y,f) 


2(d + 2£-2) 


UijL Ejj[£-|_2] (y, t) 


c\d + £ - l)(d + £ - 2)(d + 2£ + 2) 



(3.50) 


in which we denote the relevant inhnite PN series of the sonrce terms [following onr earlier 
notation (|3.30p ] by 



(3.51) 
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(3.52) 


The numerical coefficients are given by Eq. ()2.9bfl . or more explicitly 


{2k)U{2i + d){2i + d + 2) ■ ■ ■ {2i + d + 2k - 2) ' 


Notice that with our conventions the Newtonian limit, when c —> +cx), of the above defined 
relativistic moment II takes the standard Newtonian expression in any dimension d, i.e., it 
does not contain any d-dependent factors in this limit: 


4 = y d^ypVL + 0{c ^) , 


(3.53) 


where the “Newtonian” density of the fluid is p = T^^/c^. This is clear from the fact that 
the factor 2 (d- 2 ) = which appears in front of the expression of the multipole moment 
()3.50|) . cancels out precisely the d-dependent factor in the Newtonian approximation for S, 
Eq. ()3.49a|l . which is given hj T, = f p + O (c“^). 

Finally, note that the last term in ()3.48j) or ()3.5()j) is proportional to both B and £ = d —3. 
To show that this term does not contribute to II, we can Erst decompose the integral over 
d'^y in two parts: (i) an integral Xi over a compact domain r < TZ containing the two 
particles, plus (ii) an integral I 2 over the outer domain r > TZ. Even if the integration near 
the particles introduces some UV poles oc l/e,Xi will be at worst proportional to eB/e = S, 
and will [by the dehnition of the FP process, Eq. ()3.17|l ] give a vanishing hnite part at I? = 0. 
Concerning X 2 , we shall prove in Sec. IVBl below that, even if it contains infrared (IR)-type 
poles, its value is a continuous function of d. Now, because of the factor (d — 3), the value 
in 3 dimensions is zero, hence this term does not contribute to the moments and can be 
ignored in the present work (this term was neglected in Ref. Q). 


IV. SOURCE TERMS FOR THE 3PN MOMENTS IN d DIMENSIONS 


In our 3PN calculations of the gravitational wave held, we will need the expressions of 
the sources S, Sj and dehned in Eqs. above, up to orders 1/c®, 1 /c^ and 1/c^ 

respectively. The quickest method to obtain them consists in using the results of our previous 
work Q, in which the 3PN metric was expanded in terms of nine retarded potentials, 
introduced in 211 when d = 3 and generalized to d dimensions in Starting from the 
matter source densities 


a = 


(d - 2)T0° + T** 




d — 1 

we hrst dehned the “linear” potentials 


cr,: = 


<Ti, = r\ 


U = □ 


-1 

Ret 


-dvrUcr) , 


R* = □ 


-1 

Ret 


-dvrUcr,- 


(4,1) 


(4.2) 


These linear potentials were then used to construct higher “non-linear” potentials, such as 

<^kk \ 1/d — 1 


= DRet 


-47rG I aij 




d-2 


d-2 


diVd.V 


(4.3) 


and six other ones (denoted K, Ri, X, R, T) whose held equations are explicitly given 
in Eqs. (2.12) of Ref. |4^. We computed all of them for a binary system of point masses, in 
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spatial dimension d = 3 +e close to 3, at any field point in the case of linear potentials such 
as (lOl . and, for the more difficult non-linear ones like (B, in the vicinity of the particles 
as Laurent-type expansions in powers of the radial distances to them. The retardations in 
these potentials were also systematically expanded to the required PN order. 

For our present calculation of the 3PN gravitational wave field, only the expressions of 
the first seven potentials (Id, Vi, K, Wtj, Ri, X, Zij) will actually be necessary. From 
Eqs. (IXTUD and above, the sources S, Sj and may be expressed in terms of the 

“gothic” metric as 


S 



SttG 

^3 


□ 


{d - 2)h°° + R 


d — 1 


IQttG 

IbvrG 


Uh 


Oi 


nh^\ 


(4.4a) 

(4.4b) 

(4.4c) 


where G denotes by definition the gravitational constant entering the {d + l)-dimensional 
Einstein equations dSB. As underlined in Ref. Q , it is related to the usual Newton 
constant (in 3 spatial dimensions) Gw by 

G = Gw<-^ (4.5) 


where £o is an arbitrary length scale, which will enter our dimensionally regularized calcu¬ 
lation below but will drop out of the final physical observables. 

To identify the sources (lOI), it thus suffices to write Einstein’s equations R^i, = 
{SnG/— T^g^y/{d — 1)] in harmonic gauge and in terms of the gothic metric 
A possible method would be to use the expression of the Ricci tensor in terms of that we 
gave in Eq. (A9) of Ref. li^, for any dimension d. It is however quicker to use directly the 
full 3PN form of g^i, that we obtained in this reference, which can be translated in terms of 
thanks to Eqs. (A3) and (AS) of 4^ . The result not only depends on the nine introduced 
potentials [V, Vi, K, Wij, Ri, X, Zij, F), T), but the 1/c® order in g°° actually depends 
also on the 4PN (1/c®) contribution to the spatial metric gij, that was not computed in 
However, the combination entering Eq. (lOall above precisely cancels this uncomputed 
contribution, and one gets straightforwardly 


{d - 2)h00 + h" 


4r-i 


d-l 

d-2 


V^ 


8 X + 4VW + -f ^^ 
3\d-2 


d-3 

d-2 


v-^ 


K 


d-3 
d — 1 


_ sid- m -3) ; 

(d-2)2 


32T + 16 

WiMi, 


d — 1 
d-2 
4 


EX + 16VZ + 4 


d — 1 


d-2 


V^W 


2 /d-l 


'd 


+ , 

d-l 3[d-2 


V^ + 


-32 


3 \-,, 8(d-l)2(d-3) ^^^^2 


d-l 


RiVi - 


{d-2Y 


KV^ - 


8{d-l){d-3f 
{d-2f 
d-3' 




d-2 


KW 
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-16 


d-3 

d-2 


= 




VVM 


8Ri + 4 


+ 0 

d-1 

d-2 



10 




16Yi - SWi.V. + + 8 ( i —^ 1 VRi 


d-2 


+4 


d — 1 
d-2 


y. -Hd-l^KY 


{d-2y 


+ 0 



= 


1 - 


ij I 


16 


~ 2^ij^ I + ^ 



(4.6a) 


(4.6b) 

(4.6c) 


where W = Wa and Z = Zu denote the traces of the corresponding potentials. Equations 
(2.12) of Ref. 4^ then allow us to compute the d’Alembertian of these metric coefficients in 
terms of the first seven potentials, and one gets the following explicit form for the sources: 




c2 \d-2 


AttGc^ \d — 2 


d — 1 


A (1/2) 


1 


TiGd^ 


47rG a a 

:aiiV + SttG ( --) (7^1/) + dvrG 


d-2 


d — 1 


d-2) V 2 


-2V,mV - 


A{d-2y 


- iA (vw) + m.) + hzhfcl) A iKV) 

167rG 5 — d 2 8^^ jir a f~ 

■aViVi + AnG— — —CTuV^ + ^I ^—v 1 crX 


ttGc® 


d-2 


■Rg (^ 

3 \d-2 


aV^ - StiG 


{d-2y 
d-3 
d-2 


d — 1 


d-2 


aVK + —UcTiViV 


{d-l){d-2) 


+167rG ( -,—^ I aiRi — SttG— —rr^cTiiit' 


d — 1 


id-2y 


+-Wd^V + Wd^w - {dtVf - ^—^VidtVdiV 

2 ‘ 2 * 2 (d-2)3 ^ ^ {d-2f 


d-l 

d-2 


VVA^^V - 4 


d — 1 
d-2 


V,djV,d,V + 


d — 1 


(dtVi)^ + dtVdtW 


+Ad,VAW^j - AZ,A,V - + 8dy^d,R, - 2 






d-3 
d — 1 


mV, 


(4 d)(d l)(d VAdiK + 2 (^ ) W,jd,jK 


(d-2)= 


d — 1 
12 V d- 2 
1 


A(V) 


1 /d-l 


d-2 


A V^W 


2(d- 1) 


d-2^ 

Aim 


^_^-A(lE,,lE,,)-2 


V) ^ H - // 
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+4 A (r,V^ + 3) A (KV^) + (^]a ( Kw] 


+2 


S, = 


d — 1 
d-3 
d-2 
1 


A (VViVi) 


{d-2f 




d-2 




(4.7a) 


O',' 


ttGc^ 
1 /d- 1 


^ ^TrGa.y - -.— InGaVi 


d-2 


d-2 


2 \d-2 


dkVd^Vk + ll dtVd.V -\( ^ ) A (VVi) 


TiGd^ 


1 f5-d 


2 \d-2 


8{d-2y 

nGaiV^ - 2 


4 \d-2 


d-1 

d-2 


TiGaRi + 


2'kG 

d-2 


^kk^i 


+2TTGakW,k + 2nGaikVk + 2 ^^ , ^K GaVVi - 2 ^^ , ^^ Ga^K 


{d-2y 


{d-2y 


^ mV -\( mV, + ]-(^ ) dtVdtVi - 2VkdkdtV, 


4:\d-2 ' ^ 

ld{d-lf 


A\d-2 


2 V d-2 


8 (d-2)3 
d — 1 


VdtVdiV --( VidkVdkV + — ^^Vkd.VdkV 

A\d-2) 4 (d-2)2 * 


d-2 
(d- l)(d-3) 

id-2y 

1 /d- 1 

8 Id 


dkVmk - WkidkiV, + U ^ 1 dtW^kdkV - dmidkVi + dkWadiVk 


2 \d-2 

duKdiVk - (dtvm + dy/dtK) 


4(d-2)3 


l(d-l)(d-3) 


WK 


+^^(WyV,)+- 


-A (V)V 




(4.7b) 


O’. 


d-1 


' SttG \d-2 
1 r attg 


d,VdjV - - 6 ,,dkVdkV 


nGc^ \d — 
(d- l)(d-3) 


O'iV + 


d-1 

d-2 


dtV^,d,^V - d.VkdjVk - dkVAV, + 2dkV(,d,)Vk 


2(d-2)2 
1 /d- 1 


a(ya,')ir + d 




^dkVmdkV^ ^dkVmdmVk {dtVf 


2 V d-2 




+ o 


(4.7c) 


Note that although we did use the full 3PN expression of the metric in the intermediate 
steps of this calculation, the sources S, Sj and S,j actually depend only on the 2PN metric 
and on the potential Z (entering the trace of the 3PN spatial metric gij). The mass-type 
moment can now be obtained by inserting the above expressions into Eqs. (ix^-dxnn . 
thereby generalizing to d dimensions the 3-dimensional results (3.4)-(3.6) of Ref. (28| |. 

The above method to derive Eqs. dUZD not only avoids redoing some of the calculations 
of Ref. ji^, but it also yields the results in a useful form. Indeed, the Laplacians of product 
of potentials, say A(Ai?), are easier to compute than their expanded form BAA + AAB + 
2dkAdkB (where AA and AB may be replaced by their corresponding sources). In particular. 
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when computing the contributions of such Laplacians to the moment II, Eq. ()3.50|) . their 
lowest-order terms {k = 0) in Eq. (I3.51|l do not contribute to the difference between the 
dimensional and pure-Hadamard-Schwartz (pHS) regularizations; see Eq. (4.23) of Ref. [ 2 ^ 
and Sec. EB below. However, the retardation corrections {k > 1) entering Eq. (I3.51|l do 
contribute to this difference. 

To ease the reading, we classihed the various terms of Eqs. (gZTD in different sets, at each 
successive PN order: Erst the compact-support terms (proportional to a, ai or cxjj), which 
do not contribute to the difference between the dimensional and pHS regularizations below; 
second the main non-compact contributions, which are crucial for this difference; and hnally 
the non-compact terms proportional to the Laplacian of a product of potentials, which do 
not contribute to the difference at lowest order. In each set of terms, we also gathered at the 
end those which are proportional to (d — 3). These terms are absent in d = 3, and notably 
all those which involve the potential K. Note hnally that in expression (gTcl) . none of the 
terms proportional to 5ij contributes to our present calculation, since Sjj is multiplied by 
the trace-free tensor yiji in Eq. (I3.5()|) . We nevertheless quote these terms for completeness, 
as they may be useful for future works. 


V. DIFFERENCE BETWEEN DIMENSIONAL AND PHS REGULARIZATIONS 


Let us hrst recall that the general strategy we are following, in order to obtain the 
complete 3PN wave generation from two point masses, consists of two main steps. They 
have been devised at the occasion of the application of dimensional regularization (DR) to 
the problem of the 3PN equations of motion l4nl|. and are: 


(i) To obtain the expression of the 3PN mass quadrupole moment in the case of two point 
masses, using for the required self-held regularization the so-called pure-Hadamard- 
Schwartz (pHS) regularization; 

(ii) To add to the pHS result the difference between DR and the pHS regularization, which, 
as we shall see, is exclusively due to the presence of poles in d dimensions (proportional 
to l/e). 


Step (i) has already been achieved in our previous 


ja pers devoted to Hadamard- 
a |2^; the present paper deals 


regularization computations of the multipole moments 
with step (ii) of this general method, and constitutes the central part of our application of 
dimensional regularization in the problem. We refer to for a precise dehnition of the 
pHS regularization, and to for a summary and discussion of the overall method. Note 
that, in order to apply step (ii) we transformed a few terms in the expression obtained by 
inserting the ehective sources (gUl) into the multipole moments Jl so as to exactly parallel 
the form used in [ 2 ^. This is notably the case for terms that will be discussed in Sec. IVIII 
below. 

A well known result (see Refs. E M, 


21L l22j |i is that at the 3PN order, Hadamard’s 


regularization, and in fact any of its variants like the pHS one, permits the computation of 
most of the terms (both in the equations of motion and in the radiation held at inhnity), 
except for a few terms which are “ambiguous” in the sense that this particular regularization 
gives diherent results for certain divergent integrals, depending on how one performs the 
integration {e.g., by integrating by parts or not). In fact, the ambiguous integrals are those 
which exhibit some logarithmic divergencies, corresponding to the occurrence of poles in d 









dimensions. As it turns out, the structure of the ambiguous terms is always of a simple 
and limited type, and can therefore be parametrized by means of a few arbitrary unknown 
numerical constants called the ambiguity parameters'". It was shown in Refs. |^. l2^ 
that the Hadamard regularization of the 3PN mass-quadrupole moment of point particle 
binaries^^ is complete up to three and only three ambiguity parameters, which were denoted 
by K and C- 

The regularization used in the hrst work 26|| was a certain variant of the Hadamard 
regularization called “hybrid”, and the ambiguity parameters k and C, were origi nally 
dehned with respect to that hybrid regularization. The next calculation, performed in |28| . 
has been based on the pHS regularization [step (i)], and therefore we had to perform some 
numerical shifts of the values of k and C, in order to take into account the different 
reference points for their dehnition (hybrid regularization in j^, pHS one in [1^). An 
important and non trivial check of these computations has precisely been the very existence 
of a unique numerical shift for each of the ambiguity parameters, such that the results of 
both the computations 0 and [ 2 ^ are in complete agreement. Indeed, as we said, these 
two computations differ in the adopted regularizations, but they also differ by many details 
concerning their technical implementations, like the use of different “elementary” potentials. 
Indeed, in some instantaneous Poisson-like versions of the elementary potentials, say U, 
Ui, • • •, were adopted. However, in Ref. 2^ we preferred to use the retarded elementar\ 
potentials V, Vi, • • •, which are the same as in the work on the equations of motion 
and also the same as those we employ in the present paper (see Sec. HYl). 


nen tary 

Mm 


A. Difference for d-dimensional spatial integrals 

In this section we derive a general formula for the “difference” between DR and the pure- 
Hadamard-Schwartz (pHS) regularization. We shall not review the meaning and precise 
dehnition of the pHS regularization, and simply refer to Sec. HI of and Sec. IV of |2^ 
for full details. The difference investigated here concerns the typical (non-compact-support) 
terms occurring in the multipole moments, which are in the form of some spatial integrals 
over or Our investigation parallels the one of Sec. IV B in , which dealt with the 
difference for the case of Poisson and Poisson-like potentials, appropriate to the equations of 
motion. However, because the Poisson potentials depend not only on time t but also on the 
held point x, while the integrals we consider here for the multipole moments are functions 
of time t only, the derivation of the end formula will be substantially simpler than in the 
case of the equations of motion, so we shall only give the main result. 

In 3 dimensions the generic functions we have to deal with, say T’(x), are smooth on 
except at two singular points yi and y 2 , around which they admit Laurent-type expansions 
in powers (and inverse powers) of ri = |x — yi| and r 2 = |x — y 2 |.^^ When ri —> 0 we have 


The mass-quadrupole moment is the only one needed to be computed with full 3PN accuracy, thus it 
contains most of the difficult non-linear integrals, and all the ambiguities associated with Hadamard’s 
regularization. 

The function F(x) depends also on time t, through for instance their dependence on the velocities vi(t) 
and V 2 (t), but the (coordinate) t time is purely “spectator” in the regularization process, and thus will 
not be indicated. 
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(for any N E N) 


(5.1) 


F(x)= rf/,(ni)+o(rf'). 

PO 

The Landau symbol o takes its usual meaning; the coefficients fip^ni) depend on the unit 
vector Ill = |x — yi|/ri. Since the powers p can be positive as well as negative integers, the 
expansion (EH) is singular, but there is a maximal order of divergency, po G Z. 

In d dimensions, there is an analogue of the function F, which results from the same 
detailed PN iteration process as the one leading to F but performed in d dimensions (see 
the discussion in ji^); let us call this d-dimensional function F*''^)(x), where x G When 
ri —0 this function admits a singular expansion which is more complicated than in 3 
dimensions, and reads 


fplM +o{ri) , (5.2) 

Po<p<N ^ 

<io<q<qi 

with dimension-dependent coefficients /^p^g(ni) (recall that e = d — 3), and where p and q 

are relative integers whose values are limited by some po, Qq and qi as indicated. We will be 

interested here in integrands F*^'^)(x) which have no poles as £ —^ 0 (the poles in Ii being 
generated by integrating these integrands), since this will always be the case at 3PN order. 
Therefore, we deduce from the fact that F^'^)(x) is continuous at d = 3, he., lim^^s F^^'^ = F, 
the constraint 

91 

f = f pM ■ (5-3) 

9=90 

In the present paper we are interested in spatial integrals f d^xF^^\x) representing 
generic terms in the multipole moments. Here, F^'^\x.) is a non-compact-support term in 
the integrand of the multipole moments, which follows from Eqs. dUIl) in Sec. nvi above. 
We do not consider the compact support terms (proportional to a, at and aij) since their 
contribution to the moments has already been computed in Ref. and they give no contri¬ 
bution in the “difference” between the dimensional and pHS regularizations. Furthermore, 
we assume in the dehnition of the function that the derivatives of the elementary poten¬ 
tials therein are taken in the ordinary, non distributional sense (we further comment below 
on how the distributional parts of the derivatives have been taken into account into the 
formalism). Furthermore we do not need to consider here the non-compact-support terms in 
the multipole moments which have a form such that their spatial integral depends solel y on 
the boundary at infinity, |x| —>• -fcx). These terms have been discussed in Sec. IVD of j28| : 
the y p rovide a crucial contribution to the multipole moments in 3 dimensions computed 
in \2(i l28j |. However, we shall show in Sec. EH below that, thanks to the d-dimensional 
generalization of the finite part process FP^ defined in Sec. IfffAl above, these terms do not 
contribute to the difference “DR - HR” we are interested in. 

Finally, we take for a generic non-compact-support term, whose integral cannot be 
expressed as an integral at infinity, he., not of the form which is discussed in Sec. IV D of |2^ . 
The general structure of such F*-'^) is that of a multipolar factor xl times some multilinear 
functional, say V, of the elementary potentials (in d dimensions) and their derivatives, 

F('')(x) = xlV[V, V, W,j, ... , aw, • • • (5.4) 
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For the present calculation, the derivatives of potentials in this dehnition, diV,---, are 
ordinary derivatives. Many terms of Eqs. dUIl) are made of a spatial integral applied to 
some partial time derivative of a function of the type For these terms we always put 
the time derivatives outside the integral and perform Erst the spatial integral using the 
regularization, and only then apply the (total) time derivative. 

Since we shall prove in Sec. IV BI that the difference between the integrals involving 
and F^^'^ does not involve any contribution coming from divergencies “at infinity”, we limit 
ourselves to spatial integrals which extend over a finite volume in the d-dimensional space, 
say the spherical ball B(Jl) defined by |x| < TZ, where 71 denotes some arbitrary constant 
radius. The results we shall derive below will not depend on 71. In Hadamard’s regulariza¬ 
tion, and particularly in the pHS variant of it, the 3-dimensional spatial integral is defined 
by the so-called partie finie (Pf) prescription, depending on two arbitrary constants Si and 
S 2 , say 

H = Pf [ d3xF(x). (5.5) 

Of course 77 is in fact a function of time but we do not need to indicate this. By definition, 
Hadamard’s partie finie integral is given by the following limit when the radius s of two 
“regularizing volumes” surrounding the singularities tends to zero, say 

77 = lim| j d^xF(x) 

(/-3> + 1^2|. (5.6) 

p=P0 ^ ^ ^ ^ 

The symbol 1^2 means the same terms but with the singularities’ labels 1 and 2 exchanged. 
The first term represents an ordinary integral extending over the region obtained from B{7V) 
by excising two spherical balls Bi{s) and B 2 {s) centered on the two singularities, each having 
the same radius s (evidently we can always assume s -C 7Z).^^ The extra terms in (15.6j) . 
which are such that they cancel out the singular part of the first term when s —>■ 0 (so that 
the partie finie exists by definition), involve the usual (two-dimensional) spherical average 

(/>./^/(„,). (5.7) 

where dQ{ni) is the solid angle element around rii. The length scales si and S 2 (one for 
each particle) are introduced in Eq. ()5.6j) in order to adimensionalize the radius s in the 
logarithmic terms. They play a key role at 3PN order, since their appearance signals the 
presence of logarithmic divergences which correspond to poles oc 1/e in d dimensions. A 
way to interpret these constants is to say that they reflect an arbitrariness in the original 
choice of the two regularizing volumes Bi{s) and B 2 (s). 

In dimensional regularization the situation is much simpler, since the integral will be 
(so to speak) “automatically” regularized by means of the analytic continuation of the d- 
dimensional volume element. Thus, we simply have 

= [ d'^xF(‘^)(x), (5.8) 


Two balls with different radii could be used as well, without changing the results. 
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where B^'^\7V) is the d-dimensional ball with radius TZ. Given the results of the two regu¬ 
larizations, (Eini) and (lEHl), we consider what we call the difference, which is what we shall 
have to add to the pHS result in order to obtain the DR result, namely 

= (5.9) 

We shall compute VH in the limit where £ —0, keeping the pole part oc (at 3PN 
order only simple poles will occur) and the hnite term oc but neglecting 0{e). Using 
the same method as in j^, , VH can be obtained by splitting the d-dimensional integral 

(EHl) into three volumes, two spherical balls and B 2 ^\s) of radius s, which are the 

d-dimensional analogues of Bi{s) and B 2 {s), and the complementary volume in B^'^\lZ), say 
\B[^\s) UB^\s). It is clear that the integral over the latter complementary volume 
reduces, when e —0 (with hxed s), to the integral over B(7l) \ Bi{s) U B 2 {s), which is the 
hrst term in the dehnition ()5.6|) of Hadamard’s partie hnie, and does not contribute to the 
difference modulo some negligible terms 0(e). It remains thus the “local” contributions of 
the two volumes B^\s) and B^\s), which can be straightforwardly computed from inserting 
into them the local singular expansions given by Q and 1 <-»• 2. We can then connect the 
result to the corresponding result in 3 dimensions by using the constraint (jE3). Finally, we 
obtain for the difference T>H the following outcome: 


VH = 


Q 


2 +£’ 


Qi r 




1=Q0 


q + l 


-|- e Insi 


(/ + 1 " 2 + 0{e), 


(5.10) 


where the spherical average performed in d dimensions is dehned by 

(iDd i(ni) 


(/Ui - 


ft 


/(ni) ■ 


( 6 . 11 ) 


d-1 


The volume of the {d — l)-dimensional sphere, embedded into d-dimensional space, is given 
by = 27 r 2 /r(^); for instance, D 2 = dvr. Actually, we can see that the Drf_i’s cancel out 
between (ICT1D and (inii). 

Let us now comment on the inclusion in the present formalism of derivatives in a dis¬ 
tributional sense. An important feature of the pHS regularization is the systematic use of 
distributional derivatives d la Schwartz j^. It has been shown both in the contexts of the 
equations of motion HI, Hi and of the radiation held jli Hi that the purely distributional 
parts of derivatives yield a crucial physical contribution to the results at the 3PN order. In 
Hadamard’s regularization, various prescriptions are possible for the distributional deriva¬ 
tives. For instance, some generalized distributional derivatives, dehned in the extended 
Hadamard regularization j^, were used for the 3PN equations of motion in 2^. Us¬ 
ing diherent prescriptions yields diherent results, which however diher at the 3PN order by 
some terms having the form of the ambiguous terms, and therefore which merely change the 
values of the ambiguity parameters (^, n and ( in the radiation held). Now we showed 
that in DR the correct prescription for the derivatives is the one of the standard distribution 
theory This is why we have included Schwartz derivatives in the dehnition of the pHS 
regularization, which constitutes in some sense the “core” part of DR, by which we mean 
the part which computes all the difficult non-linear integrals, but leaves unspecihed a few 
terms corresponding exclusively to the “ambiguous” logarithmic divergences. Of course, 
since diherent variants of Hadamard’s regularization diher precisely in diherent dehnitions 


32 











for the ambiguity parameters, all of them could be regarded as the “core” of DR. However, 
the point is that the pHS regularization is the only one for which the final result of DR is to 
be obtained by adding exactly the “difference” in the way we have computed it in Eq. (ICTID . 
To summarize, Eq. (ICTU) as its stands is simply to be added to the pHS result, since the 
latter already includes the distributional derivatives d la Schwartz, whose contributions have 
been computed in Ref. 


B. Proof that the outer-near-zone divergencies do not contribute to the difference 

Let us recall the logic that led us to introducing and using the specihc, d-dimensionally 
generalized, hnite part (FP) process. Initially, in the MPM construction of the multipole 
expansion of the external metric, when iteratively solving Einstein’s equations, we were faced 
with some integrands Nn that had a singular behavior at the origin of the spatial coordinates, 
i.e., as r —>■ 0. One then defined the FP of the retarded integral of Nn by hrst introducing 
a factor = (r/ro)^ in the integrand, and then subtracting the “quasi-multiple” shifted 
poles — qie)~^ ■ ■ ■ (B — qk£)~^ [first term on the RHS of Eq. (jd.lbj) ]. before taking 

the continuation down to R = 0. At this stage, the integrand was, in principle, defined as a 
post-Minkowskian expansion, with good convergence properties at r —> cx), so that the poles 
(X {B — qie)~^ ■ ■ ■ {B — qk£)~^ came only from the region where r —>■ 0. Later, the external 
MPM construction was combined with a straightforward post-Newtonian (PN) iteration of 
Einstein’s equations, which took into account the interior region containing the material 
source With a generalization of the argument used in j^, one could formally relate 
the source multipole moments, used in the MPM formalism to parametrize the source, to 
integrals over the PN expansion of the effective stress-energy pseudo tensor = \g \ 
non-linear terms. 

This led to what was the starting point of our investigation, namely to formal expressions 
for the source multipole moments (of the mass type) of the symbolic form 

Jl = FP j d'^^r^XL |^+A(h)} , (5.12) 

where we recall that the overline denotes a PN (or near zone) expansion. Note that the 
presence of this PN expansion process in Eq. is crucial to its validity. Indeed, the 

argument used to derive was based on transforming MPM-expanded integrands sin¬ 

gular when r —^ 0 into PN-expanded ones diverging when r —>■ cx). (As discussed in the 
formal limit r —> cx), taken within a PN-expanded integrand, physically corresponds to the 
“outer near zone” a -C r -C A, and should not be confused with a far zone expansion r 3> A, 
in the sense of spatial infinity X°.) Technically, the transformation between the two types of 
singular integrals was based on the analytic continuation with respect to B of the integrals, 
using the fact that dr = 0. In this reshuffling from the UV (r —0) to the IR 

(r — >■ cx)) it was essential to keep the same meaning of the symbol FP in front of (EH- 
Indeed, one sees easily, by separating dr into dr and dr 

that the MPM-poles x {B — qe)~^ generated near r = 0 (when p = 1) become transformed 
in the same (modulo a sign) poles, generated near r = cxo by the singular behavior of the 
PN integrand. In addition to the poles x {B — qis)~^ ■ ■ ■ [B — qk£)~^ present in the inte¬ 
gral (ina), and generated by the behavior of the PN-expanded non-linear terms A(h) at 
r —> CX), there are also poles x associated to the singular behavior of A(h) near x = yi 
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and X = y 2 . But clearly, if we split the integral f d‘^x in a part r < TZ enclosing the two 
mass points, and a complementary part r > TZ, the latter integral will have no singularities 
associated to x = yi or x = y 2 , and therefore will have no genuine poles 

The conclusion is that the restriction of the integral to the outer near zone r > TZ 

(corresponding to the IR), say 

f d^xr^XL 

J r>Tl 

is a meromorphic function of the complex variables B and e which will have, when B and 
£ are both near zero, the same type of quasi-multiple shifted poles structure as the MPM 
quantity Fn\B) of Eq. dd.lbj] . say 


hr 


+ 


(5.13) 






{B - qie){B - q2e) ■■■{B - qks) 


+ + Cf^B + O (R2) . 


(5.14) 


The important point is that, when the expansion is written in the form the various 

coefficients C^fl, etc., are regular functions of d, which are continuous at d = 3. 

The structure proves the result we wanted, namely the fact that the IR parts of 

the integrals in d dimensions. 


= FP [/^^(5,e)] , 

B 


(5.15) 


admits when £ —> 0 the same value as the one given in 3 dimensions by the original dehnition 
of the hnite part FP (in 3 dimensions). That is to say, the two operations of taking the 
FP and the limit £ —0 commute. Indeed, the dehnition (ITT7I) of the d-dimensional FP 
operation yields 

(5.16) 

so that 

hm/^^(e:) = . (5.17) 

£—►0 

On the other hand, if we interchange the order of the two operations, we must hrst consider 
the limit when e —> 0 of (EHD, namely 

^( 3 ) 

lim/”(B,£) = ^ -^ + C® + Cf'B + O (B^), (5.18) 

and then, by applying the usual FP operation in 3 dimensions, we must discard the poles 
B~^ and evaluate the remainder at R = 0, thus 


FP 

B 


lim/^^ 


(B,e) 


= C. 


(3) 


(5.19) 


which is the same result as found in (I^TTI) . This shows that all the IR terms, formally 
depending on the boundary of the integral at inhnity, notably all those discussed in Sec. IV D 
of [ 2 ^, give exactly zero in the difference between DR and the pHS regularization. Their 
contribution to the multipole moments has already been taken into account in Ref. [ 2 ^ . 
This shows also that the outer-zone part of the last term in the expression of the moment, 
Eq. (I3.5()|l . which is proportional to Be, is actually zero in the present formalism since it is 
zero in 3 dimensions. 
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VI. COMPUTATION OF THE AMBIGUITY PARAMETERS 


As we have discussed in Sec. IV A1 the end result of the dimensional regularization (DR) 
is simply given by the sum of the pure-Hadamard-Schwartz (pHS) regularization and the 
difference" that we have investigated in the general analysis of Sec. 0 Now the pHS 
regularization of the 3PN mass dipole and qua drup ole moments of point particles binaries has 
already been computed in our previous work j28[, in which the end result of the Hadamard 
regularization (HR) was obtained as the sum of the pHS result and of some specihc ambiguity 
part parametrized by three ambiguity parameters. In the present section we construct the 
DR result and impose that it is physically equivalent to the HR one given in j^. As we 
shall show, this requirement will permit us to uniquely determine the ambiguity parameters. 


A. The 3PN mass-quadrupole moment 


Let us first state the end result of 0 concerning the 3PN mass quadrupole moment 


as computed with HR. We denote it by see Eqs. (5.9)-(5.10) in 2^ for its complete 

expression in the center-of-mass frame. In the present paper we shall not need the explicit 
formula for the moment (which includes many complicated coefficients), but simply its struc¬ 
ture, made of the sum of the pHS moment and some quite simple ambiguous contribution 
containing three and only three ambiguity parameters. The ambiguous part reads 


AT 
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f + k 


mx + m2\ (i j) 


mi 


V I A h j) 

Vi af + c ^^1 U 


+ 1 2 , 


( 6 . 1 ) 


where mi and m 2 are the two masses, y{, v\ and a\ are the position, coordinate velocity and 
coordinate acceleration of the particle 1 (1 2 denotes the same for the particle 2), and 

where the angular brackets surrounding indices mean the STF projection. All the quantities 
in (lOl are dehned in 3 dimensions; Gjsi is Newton’s constant, related to G in d dimensions 
by Eq. (EH). Obviously, since (EH) is already of order 3PN (c/. the factor 1/c®), the 
acceleration a\ is simply given by the usual Newtonian value (in 3 dimensions). 

The expression (EH) contains three ambiguity parameters {.^, k, (C}. These are the ones 
which would be dehned with respect the pHS regularization. However, the ambiguity param¬ 
eters were in fact dehned earlier in Ref. 01 which had adopted a diherent Hadamard-type 
regularization, called “hybrid”, instead of the pHS one.^® Accordingly, the pHS ambiguity 
parameters {(|, k, (f} diher from their hybrid counterparts in j^, which were denoted there 
{■Cl ^1 C}- The result, which constituted a powerful check of the computations of 0,0], is 


The hybrid regularization mainly differs from the pHS regularization in the way the “contact” (compact- 
support) terms are computed. Indeed, the hybrid regularization takes into account the so-called “non- 
distributivity” of Hadamard’s regularization, which is the fact that {F G)i f (F)i(G)i in general, where 
{F)i is the partie finie of a singular function F at the point Vi. (In this respect, the hybrid regularization 
is like the extended Hadamard regularization defined in |32|.l This introduces also some differences in 
the case of non-compact support integrals, between the “case-by-case” integration followed in and the 
systematic pHS regularization of these integrals adopted in |28| . 
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K , 


c + 
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(6.2a) 

(6.2b) 

(6.2c) 


In the present paper we prefer to stick to the original dehnition of the parameters k, C}, 
since these have already been used in the computation of the 3PN binary orbital phasing 
and in the discussion of the efficiency of the 3PN templates (see e.g. j^). Hence the 
hnal outcome from HR for the 3PN mass quadrupole moment of the binary (moving on a 
general, not necessarily circular, orbit) is written as 


4™* [’’t ’’t ’’o; 0 = h’ ’’t ’■o] + 


, 1,9 


(6,3) 


The pHS part, hrst term on the RHS, is free of the ambiguities k and but depends on the 
three regularization scales r(, rg and tq. First, ro is merely the scale we have introduced in 
the general MPM formalism, see ()3.8bjl in Sec. lHI Al and which then appears in the dehnition 
of the source multipole moments in Sec. IHI (fl This scale will disappear when we relate the 
asymptotic waveform to the local matter distribution for general extended sources. The 
other scales r[ and are specihc to the application to the case of systems of point particles 
and come from regularizing self-held ehects. By definition of the ambiguity parameters these 
scales are taken to be the same as the two scales that appear in the hnal expression of the 
3PN equations of motion in harmonic coordinates computed in Refs. j^.l2^.^^ They came 
from the regularization of Poisson-type integrals in the equations of motion, where they can 
be interpreted as some inhnitesimal radial distances used as cut-ohs when the held point 
tends to the singularities. It should be noted that r[ and are “unphysical”, in the sense 
that they can be arbitrarily modihed (though they can never be removed) by a coordinate 
transformation of the “bulk” metric outside the particles [i^, or, more consistently when 
we consider the renormalization which follows the regularization, by suitable shifts of the 
particles’ world-lines |i3. 

To get the DR result we must augment the pHS result [si, S 2 , To; k, C] computed 
for any choice of Hadamard regularization scales Si,S 2 entering Eq. (|5.6|1 . by the corre¬ 
sponding diherence Dfij^si, S 2 ; s, £ 0 ], which is made of the sum of all the contributions 
VH, Eq. (j5.1()|l . computed for all the individual non-compact support terms in the 3PN 
expression of the source quadrupole moment deduced from the explicit formulas given in 
Sec.HVl Hence this difference reads 

Vlijlsi, S 2 ; e, io] = ^ VH[si, S 2 ; e, io]. (6.4) 

non-compact 

terms in Iij 


The sum in the RHS runs over all the non-compact support terms excluding those which are 
in a form such that they depend only on the IR behavior of the integral; indeed these terms 
do not contribute to the difference (see Sec. EH). We recall also that in the calculation 


Actually only the ambiguity parameters ^ and k depend on this choice; see 


BEi 


for discussions. 
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of the difference we do not have to take into account the compact-support terms, nor the 
distributional parts of the derivatives since they are already included in the pHS result. 
In (jb.dj) we indicated that the difference depends both on the constants e = d — 3 and io 
associated with DR, and on the two scales Si, S 2 which were introduced into the Hadamard 
partie hnie (|5.ti|) . The DR result is then 

4°^^ [’^o; [-^1’ '^2, To] + Vlij [si, S2; io] ■ (6.5) 

The choice of Hadamard regularization length scales Si, S 2 in Eq. (j6.5p is arbitrary because, 
as is easily checked si, S 2 cancel out between the two terms in the RHS of ()6.5|) . so that, as 
it should be, 4^^^ depends only on the DR characteristics e and io (and also on tq which 
belongs to our general multipole moment formalism, and is in fact irrelevant for the present 
discussion). Because of this independence on the choice of the scales si, S 2 , we can choose 
them to be identical to the two specihc length scales r[, entering the 3PN equations of 
motion. Therefore we can rewrite Eq. ()6.5|) as 


4°^^ [^0; io] = [r'l, ^2, To] + Vlij [r(, r^; e, io]. (6.6) 

Let us now impose the physical equivalence between the DR result (16.61) and the corre¬ 
sponding hnal HR result (EH containing the ambiguity parameters k and (. In doing 
this identideation, we must remember, from the work on the 3PN equations of motion 
that the “bare” particle positions, and entering the DR result differ from their 

Hadamard counterparts, say and entering the equations of motion of j^, by 
some (purely spatial) shifts of the world-lines, i.e., 


ybare(^) ^ 

yb-(t) = 

These shifts have been uniquely determined in Ref. and denoted there by and ^2 (see 
Eqs. (1.13) and (6.41)-(6.43) in ji^). In the present work, we denote them by r/i and 772 in 
order to avoid any confusion with the name of the ambiguity parameter f. These shifts of 
the world-lines are crucial and must be taken into account when comparing the DR and HR 
results. Let us insist that the shifts in Eqs. EH are those which ensured the equivalence 
between the DR and HR results for the equations of motion. Having made contact in li^ 
between the renormalization scales entering the two regularization schemes in the context 
of the 3PN equations of motion, we must, by consistency, employ them to compare the DR 
and HR results for the 3PN multipole moments. The names y^^f come from the fact that the 
shifts permit to renormalize the DR result for the equations of motion, in the sense that all 
the poles oc 1/e appearing in the d-dimensional equations of motion were hnally absorbed 
into the new dehnition of the world-lines. A non trivial check of our present calculations 
will be to verify that the same shifts allow one to get finite (when e —>■ 0) hnal expressions 
for all the multipole moments, when expressed in terms of y5®2- Note that the dehnition 
of the shifts corresponds to a non-minimal subtraction. This non-minimality was needed to 
connect the DR result to the two-parameter class of HR results parametrized by arbitrary 
values of the scales r[ and (see for a discussion). Hence the shifts rji and 772 depend 
on r[ and (respectively). We shall comment more on the renormalization in Sec. IVHI 
The precise expression of the shift is 


771 (r(; £, io) 


11 G'jf ml 

y c6 



ri2 r'l 


1983 

^ 1^ 


ai, 


( 6 . 8 ) 
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together with the shift of the other world-line obtained by 1 2. Here, ai denotes the 

three-dimensional Newtonian acceleration ai = — GArm2 ni2/'rf2, where ri2 = |yi — y2| and 
ni2 = (yi — y2)/''"i2 [he., ai is the same qnantity as in Eq. ()f).l|l ]. and Gn correspondingly 
denotes the three-dimensional Newtonian constant. The expression (iniHi) seems to differ 
from the one given by Eq. (1.13) in Q , bnt this is because we have used in (jb.8|) the 
Newtonian acceleration ai in 3 dimensions, while Eq. (1.13) in Ii3 has been written with 
the help of the d-dimensional analogue a^'^^ We have 


ai-^) = 


(3+e) 


2{d-2f 
d — 1 
3 
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/C G 7712 ?"l2 ni2 , 
p/2 


= l + e 


In 


ri2q 

in 


(d=3) 


+ O 


(6.9a) 

(6.9b) 


Here we used G = GNit~\ and k = r[(d - 2)/2]/7r(‘^-2)/2 = 1 - \e\nq + 0{e‘^), where 
q = Ane^ with G = 0.577 ■ ■ ■ denoting the Euler constant. 

Evidently, since the shifts are at 3PN order, the modihcation of the mass quadrupole 
moment brought about by the latter shifts (in the sense /p[y'’®'*'®] = 7jj[y’'®“] -|- Srjlij) simply 
reads 

Srjiij = 2 mi yfr]{^ - 1 - 1 ^ 2 , ( 6 . 10 ) 

where we recall the fact that the Newtonian limit of the quadrupole in any dimension d 
takes the standard expression 1,^ = miyfy{^ -|- 1 2 -|- C>(c“^), see Eq. (I3.53|l . The 

physical equivalence between the DR and HR results simply means that we require that 
the full DR quadrupole moment, computed for the “bare” particle positions entering the 
DR delta-function source, e, io] yi®"’’®, coincides (when e — >• 0, and for the 

correct, looked-for values of k, () with the HR result r^; k, (] y™, y™]- As 

said above, the particle positions y™ entering the HR result must be identihed with the 
“renormalized” DR positions y)j®“ introduced in Eqs. (j6.7j] : y™ = Re-expressing the 
DR multipole moment in terms of the particle arguments y™ = y^®”, this requirement then 
leads to equating 


7?®'h; e. 4; y’‘", y?"] + ■ 

( 6 . 11 ) 

In other words, this equivalence is between HR and the “renormalized” result from DR. We 
hnd that the poles ~ 1/e: separately present in the two terms in the brackets of (I6.11j) cancel, 
so that the physical, renormalized, DR quadrupole moment, dehned as the RHS of (iniiD, 
is finite when e —>■ 0 and given by the limit shown. 

Let us now substitute into Eq. (IHTTD the expressions of the DR and HR quadrupole 
moments, respectively given by ()6.6jl and ra above. [We henceforth assume that on 
the LHS of Eq. (16.6|) is evaluated for y^ = Since, as we have seen, both the HR and the 

DR results have been expressed in terms of their core part, given by the pHS regularization, 
we see that, when making their comparison in (ICTll . we shall be able to remove the pHS 


r(HR) 


[r[, r2, ro; C, k, C; yr“i Tr'"] = 


£^0 . 


Note that the renormalized DR quadrupole moment is numerically equal to the original, bare quadrupole 
moment e, £o', yi“®, particular, the original, bare quadrupole moment is also finite as 

e ^ 0 (when keeping fixed y™" in taking the limit). 
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part, which is common to both sides of the equation. In this way, we obtain a relation for 
the ambiguity part A/^ of the HR quadrupole moment in terms of known quantities, viz. 


Alij 


1 9 

^ + —, K, C + - 

22 ’ ’ ^ 110 , 




io)hj 


( 6 . 12 ) 


We must now insert into the concrete result of the detailed computation of the differ¬ 

ence P/p, for all the non-compact-support terms in the explicit expression of the moment 
derived in Sec. m and following the recipe provided by Eq. dSUDD- 

The computation of P/p was performed by means of computer-aided algebraic manipu¬ 
lations, using the Mathematica software. The hnal result for P/p reads [modulo the neglect 
of 0{e) terms] 


P/p[r'i, r 2 ; £, 4] = 


Gl 


N 


mf 


+ 1 


22 220 22, /ri2r'?g3/2\ 

fc- —+ t‘” 


(i j) OD {i j) 

V i « i - ^ '^1 


(6.13) 


with the notation already used in ra and (Q. We then modify the result by including 
the effect of the particular shift which is given by Eqs. (ESI) and (IFTTUll . Thanks to this shift 
we see that the dependence of on the constants r(, r^, 5 and /q is cancelled out. More 

precisely, we hnd that the RHS of Eq. exactly takes the form of a particular instance 

of the general ambiguity term dnm), namely 




2 ) 


■^o] T ^T7(rp rj; e, 


= A/ 




9451 

9240 


, 0 , 


43 1 


330 


(6.14) 


which yields the following constraint (equivalent to three independent equations) to be 
satished by the three ambiguity parameters k, 


^lij 



c + 


9 1 
110 - 


A/ 


V 


9451 

9240 


, 0 , 


43 1 
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This immediately gives the following values for the ambiguity parameters, 

9871 

^ “ ~9240 ’ 

K = 0, 

C = -l 

^ 33 ’ 


(6.15) 


(6.16a) 

(6.16b) 

(6.16c) 


which hnally provide an unambiguous determination of the 3PN radiation held of compact 
binaries by DR. As we reviewed in the Introduction, Eqs. (16.161) represent the end result of 
DR, but in fact the results for each of the parameters k and C, have also been obtained by 
means of an independent calculation. Indeed, C = —7/33 has been shown to be a consequence 
of the Poincare invariance of the formalism (we give also an alternative, d-dimensional 
derivation of this result in Sec. IVIIII belowh the value ^ k = —9871/9240 was deduced 
from the comparison between the dipole moment and the center-of-mass position within HR 
[2^ (the latter test is equivalent to the one we shall perform below with the mass dipole in 
DR), and hnally we shall be able to check that k = 0 in Sec. IVHl ^^ 


19 


It is amusing to notice that our result for ^ happens to be related to the previous one for the equation-of- 
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B. The 3PN mass-dipole moment 


The mass-dipole moment Mi is quite interesting to consider because it satisfies a conser¬ 
vation law: Mi — Pj = 0, where Pj is the total momentum, and, as such, it can be derived 
directly from the binary’s equations of motion (instead of a wave generation formalism), as 
being linked to the conserved quantity Ki = Gi — tPi associated with the boost symmetry of 
the Hamiltonian or the Lagrangian of the binary motion. Indeed, the mass dipole moment 
is in fact nothing but the center-of-mass vector Gj ~ of the system of particles. 

Now, the center of mass vector of point particles binaries is already known at 3PN order. Its 
explicit expression was derived both in ADM coordinates Q and in harmonic coordinates 
[ 23 , Eq. (4.5) there; see also its implicit derivation in harmonic coordinates in Ref. ( 2 ^. 
[When deriving a 3PN conserved quantity we neglect the 2.5PN radiation-reaction contri¬ 
bution to the equations of motion.] We thus have the possibility of an excellent verification 
of our calculations, since the end result we shall obtain for the 3PN mass dipole moment 
Mi in DR should perfectly match with the 3PN center-of-mass Gi. In our previous paper. 
Ref. [ 2 ^, we have in fact already verified that Mi = Gi within the HR scheme, in the sense 
that we required that Mi = Gi holds, and then we deduced from this requirement the value 
of a particular combination of ambiguity parameters, namely f + k = —9871/9240. In the 
present section we shall directly show that Mi = Gi in DR, without any fine tuning of 
ambiguity parameters like in HR. 

First of all, let us recall from and the discussion in that in the present formalism 
the conserved mass dipole moment Mj is given by a slightly more complicated expression 
than the non-conserved moments II, with £> 2. Namely we have Mi = G + 6Ii, where A 
is given by the same expression as for but taken for i = 1, and where 6Ii represents a 
certain correction to it, which is given, together with the similar corrections present in the 
mass M and current dipole S'*, in Eqs. (2.22) of (in 3 dimensions). In we proved that 
the correction dli gives zero in the dipole moment at 3PN order, so that Mi = li + O (c“^). 
Now, dli is in the form of integrals at infinity (c/. the factor B in front of the integrals in 
Eqs. (2.22) of 13), and we have proved in Sec. IV HI that for such integrals the results in HR 
and DR are the same. Hence we deduce that 5Ii is also zero when applying DR and that 
Mi = Ii + O (c“^) is also true in d dimensions, modulo O (e) terms. Therefore, we need only 
discuss here the DR calculation of the main part of the dipole moment, namely Jj. 

The 3PN mass dipole moment A in HR is ambiguous, but the structure of the ambiguity 
part is very simple, as it contains one and only one ambiguity parameter f), which turned 
out to be given by the particular combination fj = ^ -\- k oi the parameters f and k which 
appeared previously in Eq. (EH). See Sec. VB in |3l for details. The structure of the 


motion related ambiguity parameter A by a simple cyclic permutation of digits: Compare 


3 ^ 


9871 

30 ^ 


with 


1987 
30 ^ ■ 


Note that the equivalence between the mass dipole moment Mi and the center-of-mass vector Gi can be 
thought of as being a consequence of the equivalence principle between gravitational and inertial masses, 
rrig = rrii. Indeed, Mi ~ while Gi ~ X ruj yb (The equivalence principle is automatically 

incorporated into the present formalism, since the motion of the point particles is geodesic, see |^.l 
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ambiguity in the dipolar case is 




2 . 


(6.17) 


Using the link we have found in (inii we can then write the HR result for the dipole moment 
in terms of the combination + k of the original ambiguity parameters in [^, hence 


r(HR) 


[r'l, rj; J + k] = 4”"®' [r'l, rQ + A/, 




( 6 . 18 ) 


which is the dipolar analogue of Eq. dnsD- However, a minor difference with (Q is that the 
3PN dipole moment happens to be independent of the cut-off scale tq. As we said above, 
the value of « could be determined in by imposing that the HR result (Q is in 
agreement with the 3PN center-of-mass position given in j23|. 

Let us now investigate what happens when using DR. Like for the case of the quadrupole 
moment, the result for the dipole moment in DR is given as the sum of the pHS dipole 
and of the difference DJj, which is made out of the sum of all the contributions of the 
non-compact support terms (excluding as usual the surface terms at inhnity) present in the 
explicit formulas of Sec. IIVI say 


VIi[si, S 2 ] e, io] = ^ VH[si, 32] e, io], (6.19) 

non-compact 
terms in li 


where si, S 2 are the two HR scales in the partie-hnie integral ()5.6|) . £ and are the DR 
scales, and each of the Dip’s are computed using Eq. (I^TTUD . Hence, 


S2] +DA[si, S2; £, 4 ] 
[r[, r'^ +DA[r(, £, 4], 


= I 


(pHS) 


(6.20a) 

(6.20b) 


where, like in the case of the quadrupole moment, we have taken advantage of the fact that 
the constants Si, S 2 cancel out from the two terms in the RHS of (j6.20a|l . to rewrite the 
result in terms of the specific length scales r(, which parametrize the 3PN equations of 
motion in 0- The last step is to renormalize the DR result by absorbing the poles in a 
spatial shift of the two particles world-lines. Of course, we must use the same shift vectors 
as in Eq. (EH), and these result in the following modihcation of the dipole moment. 


5r^Ii = mi ?7i + 1 


( 6 . 21 ) 


which is indeed checked to cancel the poles oc 1/e: of the “bare” DR dipole moment, so that 
the following limit when £ —0 is hnite. 


Mi = lira 


7i f'o] + ^T7(r',r^;e,£o)7j 


( 6 . 22 ) 


This represents our hnal “renormalized” DR dipole moment. The final result for the 

dipole moment depends on the scales and r^. We recall that this dependence does not 
correspond to any physical ambiguity, since and have the character of gauge quantities. 

Finally, after having performed the detailed calculation of the “difference” (using the same 
algebraic computer programs as for the quadrupole), and having added this difference to the 
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result for the pHS part which was obtained earlier in Ref. 2^, we found that the renormalized 
DR moment Mi given by Eq. (ESI, is in complete agreement with the conserved center- 
of-mass position Gi associated with the conservative part of the 3PN equations of motion, 
namely 

Mi[r[, r'^ = Gi[r[, (6.23) 


where Gi is explicitly given by Eq. (4.5) in Ref. jl^. We view this test as an important 
verihcation of our method and our detailed calculations. 


VII. RENORMALIZATION AND DIAGRAMMATIC APPROACH 


We have given above the hnal results obtained by combining the DR computation of the 
singular (~ 1/e) contributions to the multipole moments, coming from the vicinity of the 
point masses, with the pHS results of Ref. 2^|. This way of presenting our results is in 


close correspondence with the actual calculations we did, but it has the defect of somewhat 
hiding the logical structure of our DR results. In this section, we shall go back to basic 
methodological questions and explain in more details the logic behind DR. We shall also 
show how the examination of the structure of the DR results allows one to perform several 
checks of these results. 

Let us hrst recall that Ref. js^ presented a general method for dealing with the gravita¬ 
tional interaction of two (non spinning) compact bodies, he., bodies whose radii are of the 
same order as their gravitational radii. At the time, the main motivation for considering 
this situation was the accurate relativistic description of binary pulsar systems (he., binary 
neutron stars). Today, we have the additional motivation of accurately describing not only 
the motion but also the gravitational radiation from binary black holes (as well as binary 
neutron stars, or mixed black-hole neutron-star systems). Reference did not assume 
from the start a formal “point mass” representation of the two compact bodies but used 
instead a matching approach which combined two different approximation methods: (i) an 
“external perturbation scheme”, he., an iterative, weak-held (post-Minkowskian) approxi¬ 
mation scheme valid in a domain outside two world-tubes containing the two bodies, and 
(ii) an “internal perturbation scheme” describing the small perturbations of each body by 
the far-held of its companion. A useful outcome of this matching approach was a proof that 
to a very high approximation, the internal structures of the compact bodies were effaced 
when seen in the external scheme. More precisely, jssf (Sec. 5 there) found that the internal 
structures affected the equations of motion only starting at the 5PN level, through a term 
which is of fractional order ~ k [Gm/(c^ri 2 )]^. Here fc is a dimensionless Love number 
describing the quadrupolar deformation of one of the compact bodies under the inhuence 
of the tidal held generated by its companion. This result can be simply understood from a 
well-known Newtonian argument on the inhuence on the orbital motion of the Newtonian 
quadrupole moments induced by tidal interaction between the two compact objects (see e.g. 
Sec. 1.2 in js^). Indeed, the quadrupole moments scale as Q ~ kma^/rl 2 , where a is the 
typical size of the objects, hence in the case of compact objects for which a ~ Gmjc? we 
have in fact Q ~ {km/rf 2 ){Gm/c‘^f, which gives rise to the above mentioned correction 
to the equations of motion (and orbital phase) at the 5PN order relatively to the Newto¬ 
nian acceleration. This effacement result is the rationale for describing, up to 5PN order, 
two (non spinning) compact bodies in terms of two point masses. Technically, this means 
representing the compact bodies by a “skeleton” made of two massive world-lines, i.e., by a 
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point-particle action 


Spp = - ^maC j \J-g^,y{y^) dy^ dy’^ . (7.1) 

a 

Note that the previous reasoning suggests that, starting at the 5PN level, one will need to 
augment the effective action (HU by further terms, starting with a quadrupole-type addition 
to the monopole action (HU). At the 2.5PN level. Ref. explicitly showed how to deal 
with a point-particle description of the type (HU by using Riesz analytical continuation 
method to (uniquely) regularize the divergent integrals linked to the use of point particles in 
non-linear general relativity. It was also mentioned at the time that equivalent (2.5PN) 
results could be obtained by using an analytic continuation of the space-time dimension D, 
instead of a Riesz-type analytic continuation. 

The derivation of the equations of motion at the 3PN level turned out to be techni cally 
complicated, but conceptually satisfactory. Two independent works, 00 and 0, 
succeeded in computing using Hadamard-type regularizations most of the complicated non¬ 
linear integrals appearing at 3PN order except for a few of them, which turned out to 
be ambiguous because of the appearance of logarithmic divergencies at the 3PN order. 
Then, two further independent works, and jd^, showed that dimensional regularization 
gave unique, consistent answers, for the latter divergent integrals. A satisfactory check of 
the consistency of DR was indeed that these two independent calculations gave perfectly 
consistent final answers, though they were performed in different gauges, by completely 
different methods. 

In particular, it was found [s^ that, in Arnowitt-Deser-Misner (ADM) gauge, DR led to 
finite equations of motion (no poles oc 1/e), so that the full dynamics of the system could 
be described by an effective action obtained by adding to the d-dimensional ADM-gauge- 
fixed gravitational action the usual action for point particles coupled to gravity, namely 
Eq. (17.11) above. All the quantities appearing in this ADM plus poi nt p article action have 
finite limiting values when £ —0. By contrast, it was found in that, in harmonic 
coordinates, DR led to equations of motion containing simple poles oc 1/e, but that those 
poles could be renormalized away. There are two ways of thinking of this renormalization. 
A hrst way is to add to the usual point-particle action HU a counter-term describing a 
possible (infinitesimal) shift of the world-lines (in other words a dipole term). Then 
one shows that this dipolar counter-term is exactly what is needed to absorb the 1/e poles 
and to leave a finite answer for both the equations of motion and the bulk metric {i.e., 
the metric outside the world-lines). A second (technically equivalent) way is to use only 
the usual point-particle action HU but to consider that the “bare” world-lines entering 
HU can be decomposed in the way given by Eq. (EH), as (choosing a parametrization by 
the coordinate time, ct = yf) 


y.(t) = y’nt) + n.{t ). 


(7.2) 


where is finite as £ —> 0, but where rja, though being formally “small”, namely of 3PN 
order, contains a pole part cx 1/e: which absorbs all the poles appearing in the harmonic- 
coordinates calculations. 

Summarizing, the explicit 3PN-level calculations of the equations of motion (and of the 
pol e part of the bulk metric, see Sec. VI of 0 ) have confirmed the effacement result of 
j58l |. i.e., technically, the soundness of describing two compact bodies by the simple effective 
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action (EH). However, they also showed that, at such a high non-linearity order, it is 
crucial to use a fully consistent, and gauge-invariant regularization method. Dimensional 
regularization, which was invented precisely to preserve gauge invariance js^, Is^, , is the 
method of choice to use in this respect. 


A. Diagrammatic interpretation of the poles in DR 


As a start let us explain how one might have described the results of Sec. IVII for the 
mass multipole moments in terms of held-theory diagrams. Classical diagrammatic repre¬ 
sentations of non-linear interactions in general relativity have been introduced and used in 
several works, notably in 0, 1^ . In a previous paper of this series. Ref. 0|, we have used 
diagrams to clarify the structure of the various contributions to the equations of motion 
of two point particles. Let us do the same here for the mass multipole moments given by 
()3.50|) . 


We represent the basic delta-function sources entering as two world-lines, and each 
(post-Minkowskian) propagator as a dotted line. The various post-Minkowskian po¬ 
tentials V{x), K{x), X[x), Wij{x), etc., entering the effective sources S, S*, [see 

Sec. llVj can then be represented by drawing some dotted lines which start at the “bare” 
sources a, ai, aij, join at some intermediate vertices, corresponding to the non-linear cou¬ 
plings entering the dehnition of the non-linear potentials [such as the non-compact part of 
Wij given by ()4.3j) ]. and end at the held point x. The simpler “linear potentials”, such as 
V{x) or the “compact” part of Wij{x) (he., the part generated by a^), are just represented 
by one dotted line joining a world-line to the held point x. A product of potentials entering 
the ehective sources such as di Vj(x) dj Vj(x) is represented by juxtaposing the diagrams 
of each potential. [In this simplihed diagrammatic representation we do not explicitly in¬ 
dicate the various derivative operators which enter as “vertex factors” at the common held 
point X. However, we take care of them when they are important for the convergence prop¬ 
erties of the diagram.] Finally, we can represent the inclusion of the “multipolar factors”, 
such as xl, by adding a circled cross ®. It is then understood that one integrates over the 
“crossed vertex”, he., the held point. 

Using such a representation, the mass multipole moments are given by the sum of many 
diagrams. Note hrst that, when comparing the diagrams representing the calculation of the 
3PN multipole moment to the diagrams entering the 3PN equations of motion in US , one 
hnds that the former have a less complicated structure. Indeed, Ref. 0 has shown that 
the 3PN equations of motion involve diagrams containing up to four independent source 
points (located on the world-lines), and up to five intermediate propagators (he., hve dot¬ 
ted lines): see Figs. 2-4 in By contrast, the 3PN multipole moments only involve (if we 
treat separately, as was systematically done, the terms that can be transformed into surface 
integrals at inhnity) diagrams containing up to three source points, and four propagators. 
Examining the types of singular integrals corresponding to the possible diagrams, one then 
hnds the same rule of thumb which was found to hold in 13 for the more complicated dia¬ 
grams entering the equations of motion: namely, the only dangerously^^ diverging diagrams 
are those containing (at least) three propagator lines that can simultaneously shrink to zero 


Here and in the following we focus on the terms that generate poles oc 1/e, and we refer to them as the 
“dangerous” terms. 
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a 


b 


FIG. 1: Dangerously divergent diagrams contributing to the 3PN multipole moments. The world¬ 
lines of particles 1 and 2 are represented by vertical solid lines, the propagator by dotted lines, 
the source points by bullets, and the G symbol means a multiplication by a multipolar factor, such 
as XL, together with a spatial integration f d'^x ■ ■ ■. 


size, as a subset of vertices coalesce together on one of the world-lines. But as there are, in 
the present problem, at most three source points, this means that the dangerously divergent 
diagrams are only those represented in Fig. ^ below (or their “mirror” image obtained by 
exchanging 1 <->■ 2). 

These diagrams are also characterized by the fact that they involve, as post-Minkowskian 
diagrams (he., before explicitly performing the PN expansion, or the repeated time deriva¬ 
tives, which can introduce the acceleration of the world-line), either mf or as explicit 
factor. This reasoning is conhrmed by a scrutiny of the many explicit results reported in 
for separate pieces of the multipole moments. In the presentation of Ref. [ 2 ^ (which is less 
systematic than the more recent re-calculation of 0, but more explicit) the dangerously di¬ 
vergent integrals (in d = 3) are essentially all the terms involving the objects 
or and these terms are all multiplied by mf. Indeed, these objects are integrals of 

the type f d^xr^^ip{x) or J d^xA{r^^)ip{x), which are logarithmically divergent in d = 3, 
and lead to 1/e poles in d = 3 -|- e. 

Let us exhibit the explicit form of the terms, corresponding to the diagrams shown in 
Fig. Cl which are responsible for the poles oc 1/e in the hnal result for the multipole moments. 
Let us decompose, as in jl^, the expression for II in: (i) “hrst-order scalar” part SR 
(linear in S), (ii) second-order scalar part SIR (linear in cI/S/R), (hi) hrst-order vector 
part VR (linear in c?tSj/c^), etc. One hnds that the “dangerous” contributions to R are 
contained only in SR, SIR and VR. Moreover, one hnds that the velocity-dependent terms 
that generates poles oc l/e in intermediate calculations all cancel out in the hnal result. 
We focus here for simplicity on the non-cancelled poles, which do not depend on velocities. 
Hence, 

jdanger _ gjdanger , gjjdanger , -yjdanger f7 3 ) 

L L L L V * / 

where one checks that among the many contributions generated by inserting Eq. dUZl) into 
Eq. ()3.5()|1 the only potentially dangerous ones, in the static limit vi —0, V 2 —0, come 


Such “cancelled poles” lead to ambiguities in the finite part when working in three dimensions. This is 
taken care of in our complete results where the calculation is done in d = 3 -I- e before taking the limit 
e ^ 0. 
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gjdanger 
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- y + - 1/ - 4 Zij dij V+ 2 W,, di, K 


glldanger 
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ij 


(7.4a) 




FP — / dV^ xii -|x|" Wij dij V 


2{2i + d) TT f G dt^ 


(7.4b) 


Yjdanger 


A{2i + d-2) 
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{i + d-2){2i + d)TrfGc^ dt 


1 d{d-lf 


FP^ / XiL { - VdtVdiV 


8 (rf-2)3 




4(rf-2)^ 
■WuduV,+ 


4 (rf-2)2 


rf- 2 


rf- 1 
2{d-2) 


dtW^kdkV-d,WkidkVi+ dkWudiVA. (7.4c) 


where we set / = 

Note that the expression we used for SIIj;, in our calculations has been transformed, from 
the original form which is directly arising from the source terms given in Eqs. (EZTD above, 
by operating by parts on the terms proportional to VW and V^. This was done to exactly 
parallel the calculation of Ref. (see for instance Eq. (3.4b) there) and thereby to reduce 
the problem of evaluating the DR result to a term-by-term difference between analogous 
singular integrands. As explained in Sec. EB all the “gradient terms” generated when 
operating by parts are expressible in terms of surface integrals in the outer near-zone, and 
do not contribute to the difference between DR and pHS. We have therefore suppressed most 
of these gradient terms in Eqs. (Q, except in Eq. (I7:ial) where, as an example and as a 
reminder of the presence of such terms, we have left the terms proportional to the Laplacians 
of 1/3 and VW. 

Let us explicitly show on the example of AV^ that this term, though potentially dan¬ 
gerous, does not give rise to any pole. The linear potential V is naturally decomposed 
into V = lA -|- V 2 where lA oc mi is generated by the hrst particle, and lA oc m 2 by 
the second. In agreement with Fig. ^the dangerous contributions are cubic in mi or cu¬ 
bic in m 2 . In particular, the dangerous pieces in any term containing are (lA)^ and 
(lA)^- Let us henceforth look only at the poles generated near the hrst world-line (he., 
oc mf). In dimensional regularization, it is perfectly legitimate to integrate by parts. 
This transforms the contribution FP f d'^xr^ xi AVf into FP f d'^x A{r^ x l) Vi . Using 
A(r'® xl) = A{r^~^^ n^) = B{B + 2i + d — 2) n^, we see that the result is proportional 

to B. As we shall see in detail below the remaining integral ~ f d'^xip{x) generates a 
pole oc 1/e:. The contribution linear in AlA^ yields therefore a result proportional to B/e. 
But, by the dehnition of the d-dimensional hnite part operation, one has YV^B/e) = 0, so 
that we have indeed checked the absence of pole generated by the a priori dangerous term 
oc Al/3. Similarly for the term oc AiVW) in (I7.4aj) . 

Let us consider the various remaining terms in the integrand of SR, Eq. (TTiall . We start 
with the term oc dij V, where denotes the so-called “non compact” piece of Wj, 
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ie., the one whose source is ~ diVdjV. Again, it is easily seen that the only dangerous 
part of the integrand is ~ Vi dj Vi) dij Vi and its mirror image under the exchange 

1^2. [This term is an example of the diagram in Fig. ^.] We can compute this term by 
PN-expanding both Vi and This yields a result of the form (in the static limit) 



aoA(7f+^at5if/f + C>(!.J), 


(7.5) 


where Ui = f k G ^ with k = r(i±^)/7r^2®, is the Newtonian approximation to 
Vi, where = (Py^/dt^ is the acceleration of the hrst particle, and where and Pq are 
numerical coefficients, which depend on £. By the same reasoning as used above for the 
term cx AVp one concludes that the term cxq A Uf does not generate any pole. 

Only the second term on the RHS of m generates a pole oc 1/e which survives the 
hnite part operation. By looking at the terms contained in the last bracket on the RHS of 
Eq. ()7.4aj) one hnds that the only dangerous integrands have the same form as the second 
term on the RHS of Eq. namely proportional to ai dk Uf. Let us only give one example 
of a contribution of this form coming from the last bracket in (I7.4a|) . Consider the term 
(which can be treated to leading PN order) 

tpNC q 2 y ^ y)]NC q 2 y _ 

From the identity dkV dkV = A{V‘^/2) — VAV one has [A~^{dkV dkV)]^^ = V‘^/2, so 
that the term (uni) is of the type of d/ V. As usual the only dangerous terms are those 
proportional to d/ Vi or Vp 0/ ¥ 2 - Focusing on the hrst one, and using the fact that 

d/Vi = -a\dkVi + 0{vl), (7.7) 

one ends up with an integrand (EH) proportional to Vi a\ dk Vi or, to leading approximation 
U/ a\ dk Hi, which is indeed identical to the second term in (El). Finally, we conclude that 
the dangerous terms in Eq. Eli are of the form 

gjdanger ^ pp ^ | ~ jU | | ^ 

where Psi is a numerical coefficient which sums several similar contributions: /dsi = —Po + - ■ ■ 
[we include the factor f~^ = 1 + 0{e) into these coefficients], and where the subscript “loc” 
to the integral reminds us that one can integrate on any local neighborhood of x = yi. 

Let us now consider the dangerous terms in SIR, Eq. (I7.4bfl . In this case one must pay 
a careful attention to the dependence of the coefficients on the angular momentum index 
i. Indeed, it is important to note that there was no explicit dependence on i in Eq. EHl) 
apart from the factor xl- By contrast, the coefficients entering (I7.4bf) explicitly depend on 
i. Since fl7.4bjl has an overall factor c“®, it is sufficient to use the leading PN approximations 
for Wij and V; in view of our previous (IPN-accurate) result (j7.5j) this means that we can 
use Wij dij Vi ~ aoAUf (as usual we focus on the terms oc mf). The occurrence of an 
explicit Laplacian allows us to re-express the hrst term on the RHS of (I7.4bf) by integrating 
by parts. This leads to a term proportional to (we keep only the coefficients depending on 

^FP^/AA(|Siqx|^s 4 t/y (7,9) 
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(7.10) 


Using A (|x|^+^xl) = {B + 2){B + 2i + d) |x|^xl we get a contribution of the form 


^_FpA f 

irGc® dP ' 


B 


{B + 2){B + 2i + d) 
2i + d 


xlUI 


The pole part cx 1/e of the contribution (IT.lOj) is generated by integrating in the vicinity of 
the hrst world-line. For such a local integral the IR-converging factor |x|^ has no importance 
and we can take the analytic continuation —> 0 directly in the (localized) integrand. This 
leads to the disappearance of the ^-dependence in the factor appearing in (I71U1) . As for the 
last two terms on the RHS of ()7.4bfl . one sees that the ^dependence cancels between the 
factor cx l/{2i + d) in front, and the factors cx (2£ -|- d) multiplying the integrands xl VW 
and xlV^. Finally, we conclude that the dangerous terms in SIR are of the form 


The coefficient /3sii does not depend on i (like was the case for Psi)- The repeated time 
derivative in (17111) can then be let to act on Uf only (modulo “non dangerous” terms) 
yielding, in view of Eq. ()7.7|) , 5/ Uf ~ —a\dk Uf so that 

= / (7,12) 

A similar study of the “vector” contribution VR, Eq. (I7.4c|l . yields a result of the form 

(e + l-mf+d) L ^ 

Integrating by parts the hrst term, and taking the hnite part at R = 0 is easily seen to 
give a vanishing result [because A{xiL) = 0]. The second term of ()7.1djl . with a coefficient 
denoted /dyi, is « priori more problematic. Integrating by parts does not give a vanishing 
result (because diXn cc xl does not vanish). If present this term would have a complicated 
dependence on However, the overall coefficient /dyi of this term is the sum of many 
individual contributions, and one hnds that they all cancel out to yield /dyi = 0, so that 
hnally 

Vidanger^O. (7.14) 

The result /3yi = 0 can be obtained either by explicit calculations in d dimensions, using 
notably the explicit form of namely (to leading order) 


Wff = --(d-l){d-2)U? 


n 




+ 


Sij 


(d-i)(d-4) d{d-2y 


(7.15) 


or by considering the limiting case d = 3. In this limiting case, the poles oc 1/e are 
associated to logarithmically divergent integrals. Looking at the 3-dimensional results given 
by (8.2c), (9.3j) and (9.3k) of Ref. [1^ for VR, one indeed hnds that the terms dik 
corresponding to the d = 3 limit of the second contribution in dZH, do cancel in the 
hnal result, though the y a ppear in intermediate terms: see the hrst terms on the RHS 
of Eqs. (9.3j)-(9.3k) of 2^ with coefficients -1-2/63 and —2/63 respectively (note a small 
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misprint in (9.3k) of [ 2 ^: the overall factor nil should be understood as m\). Note also that, 
in view of the general structure (I7im derived above, it is enough to check the cancellation 
of these terms for the quadrupolar case (£ = 2) to conclude that /3vi = 0. 

Summarizing our results so far, we conclude, by adding (EHD, and (I7 :tid that the 

pole part oc ml in the £-th mass multipole moment is contained in 

(7.16) 

with a hnal coefficient (3 = (3si — /dsii- By summing the various contributions one hnds 

(3 = -^ + 0{e), (7.17) 

where the hrst term on the RHS is enough to discuss the residue of the pole cx 1/e. 


B. Renormalization of poles by shifts of the world-lines 


The result (I7:TT)D - (I7T7I1 is the explicit expression of the dangerous part of the two di¬ 
agrams of Fig. m Let us now see explicitly why it is nicely renormalized away by using 
exactly the same dipole counter-term that was found necessary in 0- The pole generated 
by (I7.1t)|l can be seen, after integrating by parts the spatial gradient (5^), as coming from 
an integral of the form 


/loc = / , 


(7.18) 


' loc 


where </3(x) is a smooth function of x (at least near x = yi). Taylor-expanding </9(x) near 
X = yi one sees that the pole in (I7.18|l comes from the zero-th term </9(yi) which multiplies 
an integral proportional to 



^3+£x^-3-3. 


Vt2+e [ 
Jo 


02+.R-'7(-2^), 


(7.19) 


where we recall that f22+£ denotes the area of the 2 + e dimensional sphere. Therefore, in 
the limit £ —> 0, the integral (I7 :thii is asymptotically equivalent to 


Jioc =-G^m?(p(yi) + C> (7) . (7.20) 

e 

This means that, when e —^ 0 the integrand Uf is asymptotically equivalent (in the formal 
sense of distributions in d-dimensional space) to 


= ^‘hG^rnldix-yi)+ 0 (e^) . (7.21) 

e 

Inserting this result into (171^ one concludes that the pole part (due to the UV divergencies 
in the neighborhoods of yi and y 2 ) of the 3PN-accurate £-th mass multipole moment is 
given by 

jpoie ^ j f^d^^Laldk |- ^ ^ (5(x - yi)| + 1 ^ 2. (7.22) 
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If we compare (I7.22|l with the leading, Newtonian approximation for namely 




d'^yiXL 1^1 ^(x — yi)| + 1 2 , 


(7.23) 


we see that the pole part i rTl i can be absorbed in a dipole-like modihcation [ afc(5(x-yi)] 
of the mass density mi 5(x — yi), or equivalently in a shift of the world-line position yi. 
More precisely, if we decompose the full yi (henceforth called the “bare” yi) as in Eq. (j7.2|l . 
with being hnite as e —> 0, but with r/i designed to absorb the pole part ()7.22j] . one 
easily checks that one needs to dehne 


in order to renormalize away this pole. Note that it was crucial to have no ^-dependence 
of the coefficients in the dangerous part (HinD in order to be able to renormalize away the 
infinite sequence of multipoles by means of the ^-independent shift rji fl7.24jl . 

In addition, by inserting the numerical value (ITTtI) of the coefficient one hnds that 
the shift needed to absorb t he p oles in the inhnite sequence of multipole moments 


ne pc 

coincides with the shift obtained in by the requirement of renormalizing both the “bulk 
metric” and the equations of motion. More precisely. Ref. found that the choice of the 
shift recalled above in Eq. ()6.8|1 [and which contains ()7.24|1 as its pole part] allowed one 
not only to get a hnite (pole-less) bulk metric and hnite equations of motion, but that the 
equations of motion coincide (when, and only when, A = —1987/3080) with the harmonic- 
gauge equations of motion, parametrized by r[ and r^, and derived using HR in Refs. 21, H- 
We recall that it is necessary to introduce some length scales r[ and associated with the 
HR of logarithmically divergent integrals in harmonic gauge. 

As we have shown here the dangerous divergencies associated with the vicinity of the 
first world-line are entirely contained in the diagrams shown in Fig. ^ and, therefore, are 
proportional to mf, without any explicit dependence on the second mass m 2 . [There is only 
an implicit dependence on m 2 via the fact that the acceleration ai is proportional to m 2 . 
But, at the level of the diagrams, ai must be considered as a pure characteristic of the hrst 
world-line.] As a consequence, we see in Eq. (I7.24|l that the dipole miT]^ needed to subtract 
the poles is also proportional to mfi This simple algebraic fact immediately leads, without 
calculations, to the result that k = 0. Indeed, the dehnition of the parameter k in Ref. [ 2 ^, 
was to parametrize a conceivable a priori ambiguity, which is indeed allowed by the weak 
assumptions of in the renormalization of the logarithmic divergencies of the type (for 
the hrst particle) 


mi 


In 


= (^ + k) m\ + Km\m 2 , 


(7.25) 


where and Si are two possible choices of regularization length scales associated to the 
hrst particle, and where we have incorporated the factor m\ associated to the divergences 
linked to yi. As fi?:^ shows, the parameter k corresponds to a mixing between diagrams 
with three legs on the hrst world-line (as in Fig. ^ and diagrams having two legs on the 
hrst world-line and one on the second. Our diagrammatic study has shown that the latter 
diagrams have no dangerous divergencies, he., that they do not introduce any conceivable 
ambiguity (even if we were working directly in d = 3, using HR). Therefore we conclude 
that K = 0. 
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The work of this section has shown that the pole in the £-th mass moment II was given 
by Eq. (TT^ whose numerical coefficients contain no dependence on the value of 1. This 
proves, in particular, that the same shift (1121, or more precisely (EHD, yields finite values 
of both the quadrupole moment {I = 2) and the dipole moment Mj [I = 1).^^ As said 
above the mass dipole moment Mj coincides with the Arnowitt-Deser-Misner (ADM) dipole 
moment or center-of-mass position Gi, such that Gi — Pit is conserved, where Pi denotes 
the total ADM linear momentum. The comparison between Mi and Gi in j2^ permitted 
to fix the value of the combination .^ + k = —9871/9240 within the HR scheme, under 
the assumption that the regularization scales si and S 2 represent some unknown but fixed 
constants, related to r[ and by some dehnite equations, and in particular take the same 
values for both the computations of the quadrupole lij and the dipole Mi. This assumption 
worked well in the case of the HR computation of the multipole moments, but failed to work 
when it was tried to assume that the same scales Si and S 2 are also those which entered 
the HR computation of the equations of motion j2^. Indeed, the work on the equations 
of motion used for the relation between si and r(, for the divergences linked to the hrst 
particle, ln(r(/si) = const + Xm 2 /rni where A was later determined to have a non-zero 
value, A = —1987/3080. Such a link is clearly incompatible with ()7.25|1 and the value we 
have found for k = 0. This means that one is not a priori allowed to assume, when using 
HR, that the scales si and S 2 represent always the same scales, fixed once and for all, and 
which can be used in different bodies of calculations. In this respect the HR is not a fully 
consistent regularization scheme. However, it can nevertheless be applied if one accepts that 
its incompleteness results in the appearance of some unknown scales Si and S 2 (generally 
in front of a few terms only), which can take different values, depending on the type of 
calculation one is doing. By contrast we have proved in Sec. I VI B1 above that the same value 
of f is consistent, in DR, with the renormalized results of both lij and Mi = Gi. This result 
constitutes evidently a solid conhrmation of the value f = —9871/9240. 


C. Comments on finite-size effects in the effective action of compact bodies 


To conclude our discussion of the diagrammatic approach to the renormalization of the 
poles which appear in harmonic gauge, let us briefly comment on the recent claim 0 that 
these poles require the introduction of new terms in the effective action describing compact 
(but extended) objects, beyond Eq. (I7.1|l and the dipole term we found above, linked to the 
shift fj7.2|l . The modified effective action proposed in Ref. has the form 


‘S'pp T ‘S'finite size , 

where Spp is the standard point-particle effective action dH} and where 




finite size 




a 


dSaRi,u{,ya)KK, 


(7.26) 


(7.27) 


with = dyf/dsa- Several claims were made in Ref. (i) that the extra terms fl7.27|l are 
necessary to “encapsulate finite size properties of the sources”, (ii) that they are linked to 


23 


Recall from Sec. I VI Bi that the conserved mass dipole moment Mi reads Mi = Ii-\-SIi, where Sf represents 
a certain correction term which, however, turns out not to contribute at the 3PN order (see |28|l. 
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the same “dangerous” diagrams that were examined in Fig. 6 of jifli l and Fig. [T] above, and 
(iii) that they entail the presence of genuine ambiguities at the 3PN level which can only 
be hxed by a matching calculation. If these statements were correct, that would mean not 
only that the basic “effacement” property (modulo 5PN-level “quadrupole-type” additional 
terms to the effective action) is incorrect, but also that the recent results, and 

this work, hxing all 3PN-level ambiguity parameters by DR are flawed. 


Let us, however, indicate why we think that the claims (i), (ii) and (iii) made in Ref. 
are not correct. First, we mention that the addition of curvature-coupling terms of the 
type indicated in (TTTTIl has already been considered in Ref. and in Appendix A of 
Ref. [^, which considered hnite-size effects in tensor-scalar gravity. Indeed, when gravity 
is partly mediated by a scalar excitation, the internal characteristics of compact objects are 
much less effaced than in the pure spin-2 case. In particular, the coupling to the spherical 
inertia moment / ~ f xcr(x) x^ can introduce extra couplings of the type of the curvature 


terms in (I7?7ll (see 
it was shown in ^ 


63l |l together with several other scalar-dependent couplings. However, 
that the use of suitable held redehnitions can transform away the 


curvature couplings (17:771) into couplings explicitly involving the gradient of the scalar held. 


f dsn Ngjifi) d,, 09 du 0 ?. As such a term does not exist in the pure spin-2 case, one sees 
that Ref. proves that can be held-redehned away. Indeed, a simple way to see it 

is to recall that the hrst-order ehect of a held redehnition of the metric + £ is 

to modify the ehective action by terms proportional to the Einstein held equations, namely 
^‘S'tot = — \ Rg^'' — ’inGT^'") eh^y [to simplify, we set the 

light velocity c = 1 here and below]. Conversely, the (a priori illicit) use of the Einstein 
held equations within an action is equivalent to a suitably dehned held redehnition e hfj,y. 
Applying this general result to (17:7111 we see that the curvature coupling terms are 

equivalent to 




finite size 


dsads'^^=d^^\y^{sa) 
m, 


+ ^®“^®ay=Wa(SaX(Sa)MaM(s'a)Mai.(s'J(5^^ni/a(Sa) -2/^(4)) (7.28) 

(where D = d+1 and = dy^/dsa) modulo a field redefinition g' = of the type 




j dSa 


[cr gf,u cR Uf, Uy\ 


-2/a(Sa)) 




(7.29) 


Here Cy, c'^, Cy are linear combinations of the coefficients c/j, cy entering ()7.27j) . namely 
= —2c'j^ = 167rG(cy — 2cR)/{d — 1), Cy = 2cy = WnGcy- After using, for instance, 
the delta function in time, d(y^(sa) — ya(Sa)), to integrate over s'^ (with the conclusion that 
s'a = Sa), one easily sees that the result (ITTHll is proportional to the Sa-integral of the d- 
dimensional delta function evaluated at a vanishing separation: S^'^^{yl^{sa) —yl{sa))- In DR, 
such a pure contact term vanishes exactly, so that we have simply = 0. [As Ref. 

uses also DR, we are entitled in using DR to discuss their claims.] Therefore we conclude 
that the proposed curvature-coupling terms (17771) are equivalent to a held redehnition of the 
type (j7.29|l . However, (17.2911 is again a “contact term” in the sense that it vanishes outside 
of the world-lines, and cannot therefore ahect the external held generated by the world- 
lines that we are interested in. In conclusion, the term (17771) can be essentially completely 
held-redehned away, and has no physical import. 
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We can give another (partial) confirmation of this result by looking at the form of the 
pole that Ref. claims to be associated with the diagrams in Fig. 6 of 113 , or Fig. m here 
(he., diagrams (c) and (d) of Fig. 7 of [H^)- Transcribing the Fourier-space result (53) of 
in a:-space, and considering the combination that enters the leading term in the multipole 
moments, one finds that, according to |^, those dangerous diagrams are equivalent, when 
£ —>• 0, to an effective mass-energy distribution of the type 


yiGold.-Roth. 


7“'00 
^(3) 


n~'ii 

^(3) 


Q G‘^ml 
e 


A(5(x-yi), 


(7.30) 


where Q is a (non zero) numerical constant and A the Laplacian. 

The result flTlTO is consistent with part of our analysis above. Indeed, using Eqs. 
dZH), (I7:ni) and (USD), our analysis has shown that the dangerous terms in the cubically 
non-linear “non compact” contributions to S, Sj and are equivalent to a term in S of 
the form 


Sfiff — — 


1 ml 


e c* 


aA6(x- yi) + 5fe'5(x - yi) 


(7.31) 


The equation (ITTinil is consistent with the first term on the RHS of (USD). But, as we 
have shown above, this term has no physical implication; only the second term, involving a 
dipole coupling a^OkS, mattered. This confirms our conclusion that the claims (i), (ii) and 


111, 


of 6 


oup 


are not correct because the terms they considered have no physical relevance. 
Note also that the “finite size” effect (unoD (formally linked to a spherical inertia moment 
/d3 

xct(x) x^, as in the tensor-scalar case of jO^) is actually a 2PN-level term. If that term 
had created physical effects linked to the finite-size of the source, this would have meant 
that the 2.5PN equations of motion 5^ had missed some 2PN violation of the effacement 
properties. As a final comment let us recall that the ADM-gauge calculations of [s^ never 
exhibited any pole. In ADM-gauge all the 3PN diagrams are finite and the whole discussion 
of possible Renormalization-Group dependent quantities evaporates away. 


VIII. QUADRUPOLE MOMENT OF A BOOSTED POINT PARTICLE 

In Sec. EH we obtained unique values for the three heretofore unknown parameters f,, 
K and (, by adding to the HR calculations of the quadrupole moment of an interacting 
binary point-mass system the additional contributions 'Dlij coming from a DR treatment 
of the singularities near yi and y 2 . In Sec. IVIII we have shown that a detailed study of the 
structure of the singular diagrams represented in Fig. fallowed one to check the values of 
both K and ^ (using information about the full computation of the dipole moment in HR to 
check the latter). Here, we shall complete our checks by giving an independent calculation 
of the third parameter (. This calculation will be based on a full DR evaluation of the 
quadrupole moment of a moving isolated particle (mi ^ 0, m 2 = 0). In another paper. 
Ref. j^, we have already checked the value of ( within a purely 3-dimensional approach, 
based on the physical situation of an isolated boosted Schwarzschild (exterior) solution with 
mass mi (and still with m 2 = 0), and without use of any self-field regularization. Therefore 
our new, DR-based, computation of ( given here can also be viewed as a further check of 
the consistency of DR. 

We thus consider the limiting case of a single particle with mass mi, moving on a straight 
line. In order to be able to discuss meaningfully this limiting case, it is important not to 
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use a center-of-mass frame for the original binary system mi, m 2 . Indeed, if we start from a 
center-of-mass frame before taking the limit m 2 —0, we shall end up with a single particle 
at rest and placed at the center of the coordinate frame used to compute the multipole 
moments. To simplify the notation, we shall suppress the index 1 on the characteristics of 
the single particle that we consider. As in we gain also some simplihcation by assuming 
that the origin of the coordinate system (with respect to which the particle is moving) which 
is used to dehne the multipole moments coincides with the position of the particle at the 
time t = 0. In other words, we consider a single particle of mass m, moving on the world-line 
j/® = vH. As was already used in jd^, the limiting case mi —>■ m, m 2 0, y\ ^ 
of the mass-type quadrupole moment of a binary system Jjj(mi,m 2 ), evaluated by HR in 
30, takes the form (at 3PN order) 


/™(m,0) = myV 


9 83 507 n® 


+ 


232 44 \ ,, , 


As we see, the ^-ambiguity enters only in a term oc We shall henceforth focus 

on this term and show how DR uniquely hxes its coefficient, he., the numerical coefficient 
C in the expression 

J,^^(m, 0) = Bmy^^y^'> + C ^ ^ . (8.2) 


To evaluate the coefficient C in DR, the first step is to obtain the D-dimensional metric, 
in harmonic coordinates, generated by a boosted point particle. We shall hrst determine 
the metric generated by a point particle at rest and then apply Lorentz invariance in D 
dimensions. There are two ways of doing this. We can start from the expressions for the 
harmonically relaxed Einstein field equations (at 3PN order) explicitly given in 0 and 
solve them by iteration, when assuming a source given by a single delta function. Another 
method consists in starting from the well-known D-dimensional Schwarzschild solution, in 
Schwarzschild-Droste coordinates, and then look for the particular harmonic coordinates 
selected by the DR treatment of delta-function sources. We have used both methods and 
checked that they fully agree. Let us indicate some details of the hrst, more pedestrian, 
approach. 

In the rest frame of a single point particle, the stress-energy tensor has = m c25W(x) 
as single non-vanishing component. This yields a “scalar” source cr(x), as used in our 
formalism, see Eq. (EH), of the form 


a 




(8,3) 


with / = 2{d — 2)/{d — l), together with cTj = 0 = cxjj. The basic scalar potential V generated 
by a, DV = —AnGa, is then found to be 


V 


f~k 


Gm 

j.d-2 ’ 


(8.4) 


where k = P[((i — 2)/2]/7r('^“^)/^. The other linear potentials are easily found to vanish, 
V) = 0, iP = 0. Going then to the various non-linear potentials, one hnds, successively, 
Ri = 0, Zij = 0, 17 = 0, as well as T = 0. Note that the vanishing of all those potentials 
results both from the treatment of contact terms in DR (namely r"(5('^)(x) = 0) and from 
the special structure of Einstein’s equations (the fact that Zij and T vanish is due to the 
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special structure of some cubic non-linearities in Einstein’s equations). Finally, besides V, 
the only non-vanishing potentials are Wij and X, which are determined by solving 


AX = WijdijV. 

As in it is usefnl to introduce the combination 


y = v- 


2 /d-3 


d-2 


4 . 16 ^ 4 ^ 

K+-X + -T = V+-X, 


which simplihes the expression of the metric. Indeed, one has 


goo = -exp ( -2^ 1 -h O 


no 


9ij = exp 


2 V 


d-2 




(8.5a) 

(8.5b) 


( 8 . 6 ) 


(8.7a) 

(8.7b) 


and goi = 0. The “gothic” metric g^^ reads, besides = 0, 

^ ^ Wkkdj + O . (8.8b) 

Note that a remarkable simplihcation occnrred in the expression (I8.8bjl of the spatial gothic 
metric. Indeed, we see form (Oi that V/c^ is proportional to Gm/c^ and therefore that 
Wij/c^ oc {Gm/c^Y while X/c® oc {Gm/c^Y- The result f)8.8bfl shows that + 

0[{Gm/cY‘^] -|- 0[{Gm/cY^] -!-■■■. The point is that there are no terms oc [Gm/c^Y 
the spatial gothic metric. One can even prove, more generally, that the spatial strnctnre 
of Einstein’s equations is snch that (for a particle at rest) contains only even powers of 
Gm/c^. The only component of the gothic metric which contains odd powers of Gm/c^, and 
in particnlar {Gm/c^Y, is the time-time component g®®, Eq. (I8.8a|l . 

By explicitly solving Eq. (I8.5a|) . we hnd, in d dimensions. 


%^-T(d-l)(d-2)V^ 


n. 


ij 


+ 


_(d- l)(d-4) d{d-2Y_ ■ 


(8.9) 


X = 


Eh 


Inserting this resnit in the RHS of (I8.5bfl then allows one to solve for X, in any dimension 
d, and we hnd 

( 8 . 10 ) 


( 8 . 11 ) 


24 \d-4: . 

Then, from Eq. (IH3I) we get, still in the rest frame: 

g®® = -A, + Gn^Y 
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where A, B and C can be expressed in terms oiV/(A and admit expansions of the type 


V 

A = 1 + Oi — + 02 + 03 — + 04 — + ... , (8.12a) 

i /2 y4 

B = 1 + 62 + ^4 + ... , (8.12b) 

i /2 p^4 

C = C 2 -^ + C4 — + ... , (8.12c) 

where, as said above, B and C contain only even powers of V/c^. The d-dependent numerical 
coefficients Oi, 02 , 03 , 62 and C 2 can be read off the results (IH31) and (Ih:hi)- (Ih:tui) above. 

It is then easy to “boost” the metric () 8 . 11|1 to a moving frame. It suffices to write it as 

= -A uV + B + uV) + C , (8.13) 


where is the H-velocity of the particle, and the unit radial H-vector orthogonal to 
the world-line.As the mass m enters only through V oc Gm, we see immediately from 
(j8.13|l that, in the “laboratory frame” where the point particle is moving, the only term in 
the gothic metric which is cubic in Gm is 




1/3 


— — 0.3 


cubic “'■5 ^6 


u>^ur 


The explicit value of the coefficient 03 in (I8.14|l is found to be 

2 


03 — 


d — 1 

d -2 


1 fd-1 


6 Vd-2 


1 /d- 1 


8 Vd -2 


d- 1 


24 V d - 4 


When e = d — 3 —>■ 0, one hnds 


as = 


l-^e + 0{e^) 


(8.14) 


(8.15) 


(8.16) 


Finally, to obtain in the lab-frame, we need to re-express the rest-frame result flHTTl) 

for V in terms of lab-frame quantities. This is simply done by saying that the rest-frame 
radial distance r entering (D can be invariantly characterized as the orthogonal distance 
r± between the world-line and the held point. In any frame, r± is given by 


r± = (>W + “p”-) (x'‘ - !/'“) {x" - y "), 


(8.17) 


where is any point on the world-line {y^ does not need to be such that — y^ be 
orthogonal to u^). Finally, we get for the part of the gothic metric deviation — 

which is cubic in Gm, 


^cubic(x, t) 



f kGm 
' ± 


3 

uV, 


(8.18) 


2"^ We have = — 1, = ’q^^n^n'' 


-bl, = 0, and 
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where = 1 /a/1 — m* = xP /c = Wj and 


rl(t) = ((5ij + UiUj) (x* - i/Xt)) (a:^ - y^{t)) , (8.19) 


where y^{t) is the point on the world-line which is lab-synchronous with the held point (at 
the same time t = y^/c). 

We have focused here on the terms cubic in Gm because, as indicated in (IHIID-dHIS) , we 
are only interested in computing the coefficient C appearing in front of the cubic term of 
(IH3). We need now to use the definition of the mass quadrupole moment fij, which is given 
by Eq. dd.bOj] . where the RHS is expressed in terms of the PN expansion of 
Introducing, as in (Id.4911 above, the notation S = ;^[(d — 2)r°° -|- r“]/c^, S* = r°Vc, 
EP = we hnally obtained the following d-dimensional expressions for the cubic terms in 
these various effective sources: 


S 


cubic 


s 


i 

cubic 


s 


cubic 


3 -I- 0{e) 1 1 

TT 1 — ^^+3e ’ 

3 + 0{e) G‘^m^ 1 

TT 1 — j ’ 

3 + 0{e) G‘^w? 1 

TT 1 — j(? ’ 


(8.20a) 

(8.20b) 

(8.20c) 


where rj_ = ri so that we have the expansion 


1 

_5+3£ 


1 / 5 -I- 3e (n\v'^Y \ 

rf+®-' V 2 ? + " y ■ 


( 8 . 21 ) 


Here ri{t) = \/— y^{t)){x^ — y^{t)) is the usual, lab-instantaneous, distance between 
the held point x and the particle y(t), and n\{t) = (x* — j/*(t))/ri. We have re-installed here 
the index 1 to distinguish the radial distance to the particle, ri = |x — y| = |x — yi|, from 
the radial distance to the origin of the lab-frame coordinate system, everywhere denoted as 
r = |x|. 

When inserting the explicit expressions (Ennii in the dehnition of the quadrupole moment, 
one ends up with a sum of d-dimensional integrals whose integrands contain several types 
of factors: an overall factor Ixl-® = |x/ro|^, various multipolar factors ~ x^, together with 
various spatial derivatives of We use dtf{ri) = —v^dif{ri) to replace time derivatives 

acting on the E’s by space derivatives. By separating the quadrupole moment in several 
contributions, as is Eq. dTHUll above, one easily checks that the leading 0{c contribution 
coming from replacing —>■ is SR and gives a vanishing contribution (after taking 

the d-modi£ed finite part). Then it takes more work to check that the 0{c~^) contribution 
VR coming from the time derivative of E* also gives a vanishing contribution. One is then 
left to evaluating an integral of the type 


R oc FP /" dR|X|^ I nf vf rf ^ jxp vf dab | • (8.22) 

J (2 2(7-|-£) J 

The dependence on £ of the global factor (not displayed here) does not import for our present 
calculation. On the contrary, the relative coefficients —(S-I-Se) and l/(7-l-£) of the two terms 
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are crucial, as there will occur below a cancellation between their lowest order contributions. 
The trick to compute Eq. ()8.22j) (for a finite value of S) is to express it, after using some 
integration by parts, in terms of parametric derivatives of “Riesz integrals”. An (Euclidean) 
Riesz integral in any dimension d is the integral 

^(a, yo, yi) = yd‘'x|x-yonx-yi|'’ = A^„b|yo-yl|“+'’+^ (8.23) 


where the numerical coefficient Nab is equal to 

r r (^) r 
r (- 1 ) r (- 1 ) r 

We hnd that we can express (Ih:^ as being proportional to 








dVody^J 


o (a o b 

oyi oy^ 


R(R + 4, -3-3e; yo, yi). 


(8.24) 


(8.25) 


Here, we have introduced, as extra parameter, the position y^ of the origin used to dehne 
the multipole moments. Up to now we have simply taken y^ = 0, but one could have dehned 
from the start the multipole moments with factors of the type |x —yo|^(a; —Inserting 
all needed factors, and explicitly evaluating the derivatives appearing in ()8.25|1 . we end up 
with a hnal answer of the type = C G'^ ml c~^ vf^\ i.e., of the form expected from 

Eq. (j8.2j) . with a numerical coefficient given, after appropriate expansion, by 


C 


FP 


B{-14:e + 9B + ■ ■■) 
7e{B - 2e) 


(8.26) 


where the ellipsis denote terms of higher order in e and/or B that do not contribute. 

We explicitly exhibit the near-hnal form (I8.26|l to emphasize the subtle nature of the 
determination of C. The result is proportional to B, which will ultimately be analytically 
continued to zero, so that one might a priori believe that C will vanish when R —> 0. 
However, this is not so because C also contains the shifted pole oc {B — 2e)~^. In addition, 
when B is non zero, (Ih:^ also exhibits a pole cx e As we explained above, the MPM 
formalism (and its subsequent PN re-expansion) imposes a specihc hnite part operation FP 
to be applied to all multipole moments. It consists in first subtracting the shifted pole terms 
and then in taking the limit B ^ 0 (see Sec. nrra . For instance, in the case of a simple pole 
of the form N{B,e)/{B — 2e), one must subtract N{2e,e)/{B — 2e) before taking R —>• 0, 
which then leads to the finite part [A^(0, e) — N{2e, e)]/{—2e). Applying this to ()8.26|1 yields 
the final result 


-2e{-lAe + 18e) _ 4 
7e{-2e) “ 7’ 


(8.27) 


which is exactly the same result as found with an independent surface-integral evaluation 
. Comparing this value to the last term on the RHS of Eq. (EH), we then conclude that 
C is uniquely fixed to the value 

C = -T. (8.28) 

in full agreement with our full two-body DR results in (I6.16|) above, and with the 
regularization-free calculations of Ref. . 
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